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INTRODUCTION 

This  program  examines  innovative  approaches  and  powerful  new  technologies  to 
identify  selective  and  potent  agents  directed  to  prevent  or  relieve  the  neuroparalytic 
toxic  actions  of  botulinum  toxin  A  (BoNTA)1.  The  focus  is  on  the  ion  channel  forming 
activity  of  BoNTs  as  a  validated  target  to  screen  for  inhibitors  of  the  translocation  of  the 
light  chain  into  the  cytosol  and  therefore  to  attenuate  the  BoNT  neurotoxicity.  The  key  of 
the  program  is  based  on  our  discovery  that  the  heavy  chain  (HC)  of  BoNT  acts  as  both 
a  channel  and  a  transmembrane  chaperone  for  the  light  chain  (LC)  to  ensure  a 
translocation  competent  conformation  during  its  transit  from  the  acidic  endosome  into 
the  cytosol  -  its  site  of  action  (Koriazova,  L.  &  Montal,  M.  Nature  Struct  Biol  10,  13-18, 
2003.  Appendix  #1)2.  This  is  an  exciting  time  to  focus  on  innovative  technologies  to 
uncover  lead  compounds  that  may  represent  a  potential  new  generation  of  useful  and 
safe  antidotes  for  BoNTs. 

Two  major  aims,  representing  different  technologies,  are  pursued: 

•  Task  1 :  Use  of  synthetic  combinatorial  libraries  (SCL)3  to  discover  specific  inhibitors 
of  the  channel  activity  of  BoNTs. 

•  Task  2:  Use  of  the  VIPR™  -Voltage/Ion  Probe  Reader4,  for  high-throughput 
screening  (HTS)  for  open  channel  blockers  of  BoNTA  channel  on 
pheochromocytoma  PCI 2  cells. 

The  combination  of  these  two  strategies  converging  on  the  channel  activity  of 
BoNT  will  optimize  the  likelihood  of  selecting  individual  compounds  from  the  libraries  in 
a  practical  time  encompassing  the  duration  of  the  grant.  Evidently,  the  two  tasks  are 
intricately  connected.  However,  inhibitor  discovery  requires  the  implementation  of  a 
robust  and  reliable  HTS.  This  has  been  the  main  focus  of  the  program  during  the  first 
year  and  is  the  subject  matter  of  this  report. 

BODY 


Our  focus  has  been  the  application  of  Fluorescence  Resonance  Energy  Transfer 
(FRET)  assays  on  nerve  growth  factor  (NGF)-treated  PCI 2  cells5  to  screen  for  blockers 
of  BoNT  HC  channel  activity.  The  HTS  assay  involves  a  Vertex,  formerly  known  as 
Aurora  Biosciences,  proprietary  phospholipid  linked  coumarin  dye  (CC2DMPE)  paired 
with  voltage  sensing  oxonol  dyes,  in  addition  to  custom  detection  machinery,  the 
VipRTM-Voltage/lon  Probe  Reader4,  as  described  in  the  application.  The  basic  concept 
underlying  this  FRET  assay  is  as  follows:  a  coumarin-linked  phospholipid-CC2DMPE, 
when  added  to  cells,  adsorbs  to  the  outer  lipid  monolayer  of  the  cell  membrane  bilayer 
and  serves  as  a  FRET  donor  to  a  mobile,  lipid  soluble  oxonol  dye,  a  bis-(1,3- 
dialkylthiobarbituric  acid)  trimethine  oxonol  (DisBAC2(3),  the  FRET  acceptor.  Excitation 
of  the  donor  results  in  FRET  provided  the  oxonol  dye  is  in  proximity  to  the  donor.  When 
cells  are  depolarized,  the  anionic  oxonol  dye  responds  to  the  transmembrane  potential 
by  diffusing  away  toward  the  inner  monolayer  of  the  bilayer,  thereby  increasing  its 
separation  from  the  membrane  anchored  donor  and  leading  to  a  decrement  in  FRET. 


4 


f 


The  ratio  change  (donor/acceptor)  is  reversed  upon  repolarization,  with  oxonol 
translocation  across  the  bilayer  in  the  ms  time  range  thereby  allowing  voltage-sensitive 
FRET  to  report  fast  voltage  changes.  For  this  assay,  cells  are  plated  into  96  well  black 
walled  tissue  culture  plates  coated  with  1  %  Matrigel  (BD  Biosciences,  growth  factor 
reduced  formulation).  Using  a  “Hi  K+  depolarization”  protocol  to  be  described,  cell 
densities  of  15  %  for  2  day  NGF-treated  PCI 2  cells  were  determined  to  be  optimal  at 
final  dye  concentrations  of  10  pM  for  CC2DMPE  and  2  pM  for  DisBAC2(3).  The  cells 
are  loaded  with  CC2DMPE  in  presence  of  Pluronic  F-127,  an  inert  polymer  which 
accelerates  loading  of  the  fatty  acid  linked  coumarin  into  the  cell  membrane  in  an 
aqueous  buffer  system.  After  the  coumarin  wash,  using  an  automatic  cell  washer 
program  designed  for  a  96  well  plate,  the  cells  are  loaded  with  oxonol  and  Pluronic  F- 
127  at  pH  7.4.  BoNT  HC  is  then  added  to  designated  wells  at  a  concentration  of  at 
least  5  pg  ml'1  to  allow  for  binding  of  the  protein  to  the  cell  outer  membrane  and 
incubated  20  min  at  22-24  °C  in  the  dark.  Next,  in  order  to  emulate  the  conditions 
prevalent  at  endosomes,  the  cells  are  exposed  to  an  acidic  environment  for  5  min  by 
adding  an  equal  volume  oxonol  containing  a  concentration  of  citrate  which  brings  the 
final  pH  to  5.5.  To  maximize  FRET,  cells  are  returned  to  pH  7.4  in  oxonol  immediately 
prior  to  the  assay  on  the  VI PR. 

To  calibrate  the  assay,  we  first  proceeded  to  use  a  “Hi  K+  depolarization”  protocol 
in  presence  of  the  well-known  channel  forming  peptide  gramicidin  A6.  Accordingly,  the 
cells  are  loaded  with  oxonol  in  a  solution  that  contains  (in  mM)  NaC1 160,  KCI  4.5,  CaCI2 
2,  MgCI2  1,  Hepes  10,  Glucose  10  (pH  7.4  with  NaOH).  Depolarization  in  “Hi  K+”  buffer 
that  contains  (in  mM)  KCI  164.5,  CaCI2  2,  MgCI2 1,  Hepes  10,  Glucose  10,  (pH  7.4  with 
KOH)  results  in  an  increase  of  the  coumarin-to-oxonol  signal  ratio  as  the  oxonol  moves 
away  from  the  coumarin  and  FRET  decreases,  i.e.  coumarin  signal  increases  while 
oxonol  signal  decreases.  The  final  data  are  expressed  as  normalized  ratio  between  the 
two  dyes  in  a  column  by  column  display,  as  shown  in  the  example  displayed  in  Fig.  la 
(Control).  In  presence  of  10  pM  gramicidin,  the  signal  is  attenuated  or  nearly  abolished 

(Gramicidin),  in 
agreement  with 
expectations  that  the 
K+-conducting 
gramicidin  channel 
would  collapse  the 
membrane  potential 
generated  by  the 
activity  of 
endogenous  K+ 
channels.  This 
calibration  defines 
an  end  point  with 
which  to  compare 
the  anticipated  effect 
of  BoNT  HC 
channel. 


a  b 
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For  the  HTS-FRET  assay  of  the  HC,  we  have  exploited  the  fact  that  the 
endogenous  channels  of  PC12  cells  are  essentially  impermeant  to  Cs+  whereas  the  HC 
channel  is  highly  conductive  for  Cs+,  as  shown  in  independent  measurements  using 
patch  clamp  recordings7.  Accordingly,  the  cells  are  loaded  with  oxonol  in  a  solution  that 
contains  (in  mM)  NaC1 160,  KCI  4.5,  CaCI2  2,  MgCI2 1,  Hepes  10,  Glucose  10  (pH  7.4 
with  NaOH),  and  stimulated  in  a  Hi  Cs+  medium  that  contains  (in  mM)  CsCI  164.5, 
CaCI2  2,  MgCI2  1,  Hepes  10,  Glucose  10  (pH  7.4  with  CsOH).  Representative  results 
are  shown  in  Fig.  1  b.  An  increase  in  fluorescence  ratio  was  determined  in  presence  of 
the  BoNT  HC;  this  is  consistent  with  the  fact  that,  under  these  conditions,  the  BoNT  HC 

channel  is  the  major  channel  which  conducts  Cs+. 

The  conditions  developed  thus  far  require  optimization  which  is  a  major  currrent 
effort  in  the  laboratory;  however,  the  initial  results  validate  the  HTS-FRET  assays  on 
NGF-treated  PCI 2  cells  to  screen  for  blockers  of  BoNT  HC. 

KEY  RESEARCH  ACCOMPLISHMENTS 

•  A  key  step  in  the  intoxication  by  BoNT  is  the  translocation  of  internalized  toxin 
across  intracellular  membranes  to  reach  its  cytosolic  targets.  A  fundamental 
discovery  was  the  demonstration  that  the  heavy  chain  acts  as  both  a  channel  and  a 
transmembrane  chaperone  for  the  light  chain  protease  to  ensure  a  translocation 
competent  conformation  during  transit  from  acidic  endosomes  into  the  cytosol 
(Appendix  #1).  Thus,  the  stage  is  set  for  pursuing  the  identification  of  channel 
blockers  specific  for  the  BoNT  heavy  chain. 

•  A  high-throughput  screen  exploiting  the  VIPR™-Voltage/lon  Probe  Reader,  is  now  a 
proven  strategy  for  the  search  of  compounds  that  modulate  or  block  the  channel 
activity  of  BoNT  on  NGF-treated  PCI 2  cells. 

REPORTABLE  OUTCOMES 

•  Publications 

Appendix  #  1 :  Koriazova,  L.  and  Montal,  M.  Translocation  of  botulinum  neurotoxin  light 
chain  protease  through  the  heavy  chain  channel.  Nature  Struct.  Biol.  [Published  online 
2  December  2002;  doi:10.1038/nsb879]  10:13-18  (2003). 

Appendix  #  2:  Harms,  G.S.,  Orr,  G.,  Montal,  M.,  Thrall,,  B.D.,  Colson,  S.D.,  and  Lu, 
H.P.  Probing  conformational  changes  of  gramicidin  ion  channels  by  single-molecule 
patch-clamp  fluorescence  microscopy.  Biophys.  J.  85:1826-1838  (2003). 

•  Invention  Disclosure 

UCSD  #SD2003-007  entitled:  Method  for  antidotal  therapy  for  botulinum  toxin  (brief). 
Title:  The  botulinum  neurotoxin  channel  as  a  validated  target  for  intervention  aimed  to 
inhibit  the  translocation  of  the  light  chain  into  the  cytosol  and  therefore  to  attenuate  the 
neurotoxicity  of  botulinum  neurotoxins. 

Inventor:  Mauricio  Montal 
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CONCLUSIONS 


The  key  discovery  is  that  the  heavy  chain  of  BoNT  acts  as  both  a  channel  and  a 
transmembrane  chaperone  for  the  light  chain  to  ensure  a  translocation  competent 
conformation  during  its  transit  from  the  acidic  endosome  into  the  cytosol  -  its  site  of 
action.  The  light  chain  is  a  Zn2+-metalloprotease  that  cleaves  the  protein  components 
involved  in  synaptic  vesicle  fusion  with  the  neuronal  membrane,  thereby  abrogating 
synaptic  transmission1.  To  accomplish  this  task,  the  heavy  chain  operates  as  a 
transmembrane  protein-conducting  channel2:  the  channel  is  occluded  by  the  light  chain 
during  transit,  and  open  after  completion  of  translocation  and  release  of  cargo,  acting 
intriguingly  similar  to  the  protein-conducting/translocating  channels  of  the  endoplasmic 
reticulum  (ER),  mitochondria,  and  chloroplasts.  This  finding  has  outlined  a  novel  way  of 
thinking  about  BoNT  neurotoxicity  shifting  the  focus  of  attention  on  its  translocation 
within  cells  rather  than  on  the  protease  activity  of  the  light  chain,  which  is  known  not  to 
be  toxic  unless  it  is  internalized. 

Current  progress  focused  on  screens  for  blockers  of  the  BoNT  HC  channel  as  a 
target  for  lead  discovery  are  exciting  and  encouraging.  A  major  outcome  of  the 
proposed  program  would  be  a  blueprint  to  uncover  blockers  of  the  protein-conducting 
channel  of  BoNT  which  may  be  developed,  in  due  turn,  into  prophylactic  agents  to  be 
administered  to  military  personnel  prior  to  deployment  or  to  civilian  populations  at  risk. 
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Translocation  of  botulinum  neurotoxin  light 
chain  protease  through  the  heavy  chain 
channel 

Lilia  K.  Koriazova  and  Mauricio  Montal 

Published  online  2  December  2002;  doi:10.1038/nsb879 


Clostridia I  botulinum  neurotoxins  (BoNTs)  abort  the  process  of  neurotransmitter  release  at  presynaptic  motor 
nerve  terminals,  causing  muscle  paralysis.  An  enigmatic  step  in  the  intoxication  process  is  the  mechanism  by 
which  the  neurotoxin  heavy  chain  (HC)  forms  the  conduit  for  the  translocation  of  the  light  chain  (LC)  protease 
across  the  endosomal  membrane  into  the  cytosol,  its  site  of  action.  Here  we  investigate  the  mechanism  of  LC 
translocation  by  using  the  combined  detection  of  channel  currents  and  substrate  proteolysis,  the  two  hallmark 
activities  of  BoNT.  Our  data  are  consistent  with  the  translocation  of  the  LC  through  the  HC  channel  and  show  that 
the  LC  protease  activity  is  retrieved  in  the  trans  compartment  after  translocation.  We  propose  that  the  BoNT 
HC-LC  complex  embedded  in  the  membrane  is  a  transmembrane  chaperone,  a  dynamic  structural  device  that 
prevents  aggregation  and  achieves  translocation  of  the  LC.  In  this  regard,  the  complex  is  similar  to  the  protein 
conducting/translocating  channels  of  the  endoplasmic  reticulum,  mitochondria  and  chloroplasts. 


Clostridial  botulinum  neurotoxins  (BoNTs)  act  as  sequence- 
specific  endoproteases  that  cleave  unique  constituents  of  the 
SNARE  complex  —  the  synaptic  vesicle  doddng/fusion  com¬ 
plex1”6.  A  widely  held  view  is  that  BoNTs  enter  cells  via  receptor- 
mediated  endocytosis2.  Exposure  of  the  BoNT  holotoxin  to  the 
acidic  pH  of  the  endosomes7’3  induces  a  conformational  change 
of  the  disulfide-linked  dichain  toxin,  allowing  the  passage  across 
the  membrane  of  the  light  chain  (LC)  protease  through  die  puta¬ 
tive  heavy  chain  (HC)  channel  and  into  the  cytosol.  Until  now, 
there  has  been  no  direct  evidence  for  the  translocation  of  the  LC 
through  the  HC  channel.  The  crystal  structures  of  BoNT  A9  and 
BoNT  B10  confirmed  a  tri-modular  design  of  the  neurotoxin 
protein:  a  receptor  binding  domain,  a  translocation  domain  and 
a  catalytic  domain.  However,  these  structures  were  determined 
at  pH  7  and  pH  6  for  BoNT  A  and  B,  respectively,  and  represent 
a  single  conformation  of  the  inactive  holotoxin  in  the  soluble, 
but  not  the  membrane-embedded,  form.  They  also  provide  no 
compelling  molecular  details  about  the  elusive  translocation 
process. 

Here,  we  use  the  combined  detection  of  channel  currents  and 
substrate  proteolysis,  the  two  hallmark  activities  of  BoNTs,  to 
determine  if  HC  channel  activity  and  LC  translocation  are  corre¬ 
lated.  We  demonstrate  that  channel  activity  requires  acidifica¬ 
tion  of  the  medium  in  the  aqueous  compartment  to  which 
BoNT  A  is  added  (ris)  and  a  reducing  environment  in  the  oppo¬ 
site  compartment  (trans).  Under  these  conditions,  we  detected 
the  LC  protease  in  the  trans  compartment,  consistent  with  its 
translocation  through  the  HC  channel.  CD  data  indicate  that 
only  die  unfolded  conformation  of  LC  A,  the  LC  of  BoNT 
serotype  A,  correlates  with  both  channel  and  protease  activities 
of  BoNT  A.  Together,  the  results  imply  that  the  LC  partially 
unfolds  at  the  endosomal  pH,  passes  through  the  HC  channel, 
refolds  at  the  interface  and  dissociates  from  the  HC  in  the  reduc¬ 
ing  cytosolic  milieu  where  it  cleaves  its  substrate  SNAP-25.  These 


findings  suggest  a  novel  modality  of  transmembrane  chaperone 
action  based  on  protein  conducting/translocating  channels. 

LC  ocdudes  the  HC  channel  in  unreduced  BoNT  A 

To  investigate  the  mechanism  of  LC  translocation,  we  relied  on 
the  channel  activity  of  BoNT  A  holotoxin  and  HC  A  in  lipid 
bilayers  and  used  single  channel  recordings11”13  as  a  sensitive 
readout  to  monitor  whether  the  LC  blocks  the  HC  channel  as  it 
passes  from  the  cis  compartment  of  the  membrane  to  the  trans 
one.  Evidence  that  BoNT  A  forms  channels  in  lipid  bilayers16’15 
and  PC12  cells16,  predominantly  after  exposure  to  an  acidic  pH 
and  only  after  chemical  reduction,  has  led  to  the  hypothesis  that 
the  HC  forms  channels.  Our  protocol  here  was  to  add  BoNT  A  to 
the  cts  solution,  which  was  held  at  pH  7.0  to  enhance  binding  to 
the  gangliosides-containing  bilayers.  After  2  min,  the  cis  solution 
was  acidified  (4.5  <  pH  £  5.5)  while  the  trans  compartment  was 
held  at  pH  7.0.  Under  these  conditions,  unreduced  holotoxin 
BoNT  A  does  not  form  channels  (Fig.  la).  In  contrast,  unre¬ 
duced  native  (Fig.  lb)  and  recombinant  (Fig.  lc)  HC  from  BoNT 
A  forms  channels  similar  to  those  recorded  with  reduced  BoNT 
A  (Fig.  Id).  Clearly,  the  HC  is  a  channel  irrespective  of  redox 
state.  These  results  raised  the  intriguing  question  of  whether,  for 
the  unreduced  holotoxin,  the  HC  channel  is  occluded  by  the 
anchored  LC. 

Disulfide  reduction  unblocks  the  HC  channel 

To  test  this  inference,  we  examined  the  effects  of  the  membrane 
nonpermeable  reductant  tris-(2-carboxyethyl)  phosphine17 
(TCEP)  on  BoNT  A  channel  activity.  The  most  telling  results 
emerge  by  examining  the  contrasting  findings  generated  by  the 
location  of  TCEP  in  the  ris  or  the  trans  compartments  (Fig.  1). 
Unreduced  holotoxin  (data  not  shown),  or  BoNT  A  in  the 
absence  of  TCEP  in  the  acidified  ris  compartment,  does  not  form 
channels  (Fig.  la).  Presence  of  TCEP  in  the  trans  compartment. 
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however,  elicits  channel  activity  (Fig.  Id).  When  reduced  holo- 
toxin  is  added  to  the  as  compartment,  followed  by  acidification, 
channels  are  formed  and  there  is  no  need  for  TCEP  in  the  trans 
compartment  Notably,  addition  of  TCEP  to  the  cis  compart¬ 
ment  1  min  after  acidification  fails  to  evoke  channel  activity 
(Fig.  la).  The  absence  of  channel  activity  when  TCEP  is  added  to 
the  cis  compartment  after  acidification  is  interpreted  as  indica¬ 
tive  of  the  translocation  of  the  LC  through  the  HC  channel,  ren¬ 
dering  the  critical  disulfide  link  between  LC  and  HC  inaccessible 
to  the  membrane-impermeant  TCEP.  Presumably,  the  LC 
remains  attached  to  the  HC,  occluding  the  HC  channel,  whereas 
reduction  in  the  trans  compartment  releases  the  LC,  thereby 
unblocking  the  HC  channel.  This  observation  implies  that  the 
LC  blocks  the  HC  channel  during  its  translocation. 

The  LC  selectively  blocks  the  HC  channel 

This  prediction  was  tested  using  isolated  HC  and  LC  and 
addressing  the  specific  effect  of  the  LC  on  HC  channels  depend¬ 
ing  on  its  location  in  the  cis  or  the  trans  compartments  (Fig.  2). 
Neither  native  nor  recombinant  LC  A  forms  channels14’15.  A  ten¬ 
fold  molar  excess  of  LC  in  the  acidified  cis  compartment,  in  the 
absence  of  reductant,  does  not  affect  the  HC  channel  activity 
(Fig.  2a).  In  contrast,  a  stoichiometric  amount  of  LC  in  the  neu¬ 
tral  trans  compartment  attenuates  and  eventually  silences  the 


HC  channel  activity  (Fig.  2b):  the  single  channel  conductance 
(y)  is  sharply  reduced  (Fig.  2b  inset)  from  110  pS  to  25  pS,  and 
the  channel  openings  become  brief,  appearing  as  transient  spikes 
at  this  time  resolution.  This  pattern  is  characteristic  of  open 
channel  blockers  and  arises  from  the  progressive  shortening  of 
the  residence  time  in  the  open  state12.  Blocking  is  elicited  by  both 
native  and  recombinant  LC  on  both  the  native  and  recombinant 
HC  channels,  is  independent  of  the  polarity  of  the  applied  volt¬ 
age  and  is  not  produced  by  thermally  or  guanidine-HCl- 
denatured  and  inactivated  LC  (Fig.  2c).  Thus,  at  pH  7.0,  the  LC 
selectively  blocks  the  HC  channel. 

LC  translocation  revealed  by  Its  protease  activity 

To  examine  if  the  LC  protease  goes  through  the  HC  channel,  we 
developed  a  high-sensitivity  ELISA  and  scaled  up  the  single 
channel  measurements  to  detect  numerous  (>1,000)  channels 
(Fig.  3a).  The  trans  compartment  contained  the  LC  substrate, 
SNAP-25  (refs.  2,3).  A  period  of  15-30  min  was  fixed  for  each 
assay;  afterward,  the  entire  volume  of  the  trans  chamber  was  col¬ 
lected  and  assayed  for  SNAP-25  proteolysis.  No  cleavage  of 
SNAP-25  was  detected  in  the  absence  of  BoNT  A  (control)  or  in 
the  presence  of  LC  A  in  the  cis  compartment  (Fig.  3c).  The  latter 
control  indicated  that  the  bilayer  prevented  mixing  of  isolated 
LC  (cis)  with  SNAP-25  (trans)  and  confirmed  that  LC  perse  can- 
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Fig.  2  Blocking  of  HC  A  channels  by 
the  LC  A.  a,  Single  channel  currents 
of  native  BoNT  HC  in  the  presence  of 
LC  in  the  cis  chamber  LC  A  (25  pg 
mh1)  was  added  after  the  onset  of  HC 
A  channel  activity.  For  this  record,  y  = 
124  pS  and  P0  =  0.56.  The  mean  y  e 
113  ±  9  pS  and  n  =  5.  I>,  Single  chan¬ 
nel  currents  of  native  BoNT  HC  in  the 
presence  of  LC  in  the  trans  chamber. 
LC  A  (2.5  pg  mH)  was  present  from 
the  beginning  of  the  experiment 
The  segment  delimited  by  the  arrows 
is  displayed  at  higher  resolution.  For 
this  record,  y  =  29  pS  and  0.51. 
The  mean  y  =  25  ±  4  pS  and  n  =  5. 
Note  the  occurrence  of  fast  transi¬ 
tions  between  open  and  blocked 
states,  c,  Single  channel  currents  of 
native  BoNT  HC  in  the  presence  of 
denatured  LC  A  (10  pg  mM)  in  the 
trans  chamber.  LC  A  was  inactivated 
by  incubation  in  6  M  guanidinium-CI 
for  30  min  before  the  addition. 
Inactivation  was  confirmed  by  direct 
assay  of  SNAP-25  cleavage.  For  this 
record,  y  =  123  pS  and  P0=  0.54.  The 
meany=112±9  pSandn  =  5.  Protein 
concentration  was  5  pg  mh1  for  HC 
Other  conditions  were  as  in  Fig.  la. 


not  permeate  through  an  unmodified  lipid  barrier.  Qeavage  of 
SNAP-25  requires  the  presence  of  TCEP  on  the  trans  side  and 
pH  5.0  on  the  cis  side  (Fig.  3c),  a  condition  that  correlates  with 
the  insertion  of  multiple  channels  in  the  bilayer  (Fig.  3a,  equiva¬ 
lent  to  Fig.  la  for  single  channels).  This  correspondence  argues 
that  there  is  proteolytic  activity  of  the  LC 


established  a  calibration  curve  (Fig.  3d).  The  minimum  BoNT  A 
LC  activity  was  discerned  at  3  x  10-18  M,  producing  a  10  ±  5% 
cleavage  of  SNAP-25.  At  3  x  10“17  M>  25  ±  5%  cleavage  was  reli¬ 
ably  measured.  This  value  approximates  the  value  detected  with 
TCEP  on  the  trans  side  and  at  pH  5.0  on  the  cis  side  (Fig.  3c), 


on  die  trans  side  only  under  conditions 
in  which  the  channel  activity  of  the  holo- 
toxin  is  detected.  This  is  consistent  with 
the  concomitant  absence  of  channel 
(Fig.  3 b)  and  LC  protease  (Fig.  3c)  activi¬ 
ty  when  the  pH  on  the  cis  side  was  4.5. 

To  estimate  the  number  of  LC  mole¬ 


cules  translocated  across  the  bilayer,  we 


Hg.  3  Detection  of  the  endopeptidase  activity 
of  translocated  LC  A  in  bilayer  experiments, 
a.  Time  course  of  insertion  of  multiple  BoNT  A 
channels  into  a  membrane  at  pH  5.  Records 
obtained  at  -10  mV.  The  cis  solution  was  0.1  M 
Kd,  1  mM  Cad2  and  5  mM  HEPES,  pH  7.1  ± 
0.1.  After  addition  of  BoNT  A  (0.25  pg  mh1) 
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the  pH  was  adjusted  to  5.0.  The  trans  solution 
was  0.01  M  KC1,  1  mM  CaCI*  5  mM  HEPES, 
pH  7.1  ±  0.1,  4  pM  N-terminal  biotinylated 
SNAP-25,  0.5  mM  TCEP  and  10  pM 
Zn(CH3COO)2.  b,  No  detectable  channel  activ¬ 
ity  of  BoNT  A  at  pH  4.5.  Conditions  identical 
to  (a)  except  that  the  pH  in  the  ds  compart¬ 
ment  was  4.5.  c;  SNAP-25  cleavage  in  samples 
collected  from  the  trans  chamber  of  bi layer 
experiments  equivalent  to  the  representative 
records  shown  in  (a)  and  (b).  Error  bars  indi¬ 
cate  s.d.  and  n  *  6  experiments,  except  for 
BoNT  A  (+TCEP)  at  pH  5.0  (a),  for  which  n  =  10. 
Probabilities  are  indicated  by  the  asterisk 
from  t-tests  performed  with  respect  to  control 
for  each  experimental  condition  separately: 
'n.s/  denotes  not  significant;  asterisk,  P  < 
0.001.  d,  SNAP-25  cleavage  as  function  of 
BoNT  A  concentration  determined  by  EUSA 


(n  =  4).  The  values  obtained  for  the  isolated  LC 
A  at  selected  concentrations  agree  well  with 
those  shown  in  (d). 
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which  corresponds  to  -3,600  LC  molecules  in  the  experimental 
samples.  This  number  is  in  dose  agreement  with  the  maximum 
number  of  channels  inserted  in  the  bilayer  under  these  assay 
conditions,  -2,000  (Fig.  3  a),  as  estimated  from  the  macroscopic 
membrane  conductance  -1  x  10s  pS  and  y  =  50  pS  measured  in 
0. 1 M  KC1.  Together,  these  findings  demonstrate  translocation  of 
the  LC  through  the  HC  channel  and  retrieval  of  the  LC  protease 
activity  at  neutral  pH. 

How  does  the  LC  go  through  the  HC  channel? 

The  diameter  of  the  BoNT  A  or  HC  A  channel,  estimated  from  y, 
is  -15  A  (refs.  14,18),  raising  the  question  of  how  a  pore  of  this 
size  could  accommodate  the  globular  LC  with  dimensions  of 
55  A  X  55  A  x  62  A  (ref.  9).  These  considerations  imply  that  the 
LC  unravds  at  the  endosomal  pH  (5.1  <pH<  5.4)7  to  go  through 
the  channel  and  refolds  after  exiting  at  the  neutral  cytosolic  pH. 
To  assess  the  pH-induced  structural  changes  of  the  LC  and  HC 
assodated  with  their  respective  protease  (Fig.  4a)  and  channel 
(Fig.  4b)  activities,  we  used  CD  spectroscopy  (Fig.  4c).  At  pH  7.0, 
the  a-helical  content  of  the  individual  LC  A  and  HC  A  is  compa¬ 
rable  to  that  calculated  from  the  crystal  structure  of  BoNT  A9. 
The  a-helical  content  of  the  HC  A  was  unaffected  by  pH.  For 
LC  A,  there  is  a  sharp  decrease  in  helicity  at  pH  5.0  (Fig.  4c), 
which  correlates  with  the  occurrence  of  channel  (Figs.  1  a,  3a,  4b) 
and  protease  (Figs.  3c,  4a)  activities  of  BoNT  A  In  contrast,  at 
pH  4.5,  there  is  a  drastic  increase  in  LC  helicity  coincident  with 
the  absence  of  channel  (Figs.  3b,  4b)  and  protease  (Figs.  3c,  4a) 
activities.  Together,  the  data  indicate  that  only  the  unfolded  con¬ 
formation  of  LC  A  correlates  with  both  channel  and  protease 
activities  of  BoNT  A  This  suggests  that,  under  the  conditions 
prevalent  in  the  cell  after  BoNT  A  entry — specifically,  the  endo- 
some/cytosol  pH  gradient7  and  a  cytosolic  reducing  milieu19  — 
the  LC  unfolds  or  partially  unfolds  in  the  endosome,  goes 
through  the  HC  channel,  dissociates  and  refolds  at  the  interface, 
and  cleaves  SNAP-25  in  the  cytosol.  At  pH  4.5,  the  acid-induced 
conformational  change  of  BoNT  A  allows  insertion  of  the  HC 


Fig.  4  BoNT  A  endopeptidase  activity  correlates  with  BoNT  A  channel 
activity  and  the  unfolding  of  LC  A.  a,  Endopeptidase  activity  in  samples 
collected  from  the  trans chamber  of  bilayer  experiments  (Fig.  3c)  as  func¬ 
tion  of  pH  in  the  ds  compartment  The  trans  compartment  pH  was  7.0. 

b,  Channel  activity  of  BoNT  A  and  HC  A  as  function  of  pH  in  the  ds  com¬ 
partment;  the  trans  compartment  pH  was  7.0.  Absence  or  presence  of 
channel  activity  is  arbitrarily  defined  as  0  or  1;  n  =  6  for  HC  A  and  n  *  10 
for  BoNT  A  (Fig.  3a, c).  Other  conditions  equivalent  to  those  of  Fig.  3a. 

c,  a-helical  content  of  LC  A  and  HC  A  as  function  of  pH  calculated  from 
far  UV-CD  measurements  carried  out  at  25  *C  (n  =  3). 

into  the  membrane,  but  the  high  helical  content  of  the  LC  traps 
it  within  the  HC  channel  and  occludes  the  channel,  thereby  pre¬ 
cluding  detection  of  channel  or  protease  activities. 

The  HC  acts  as  a  protein-translocating  channel 

The  interplay  between  an  unfolded  LC  embedded  within  the  HC 
is  reminiscent  of  the  maintenance  of  an  unfolded  or  partially 
folded  state  of  polypeptides  by  chaperones.  Is  the  HC  acting  as  a 
transmembrane  chaperone  for  the  LC  to  ensure  a  translocation- 
competent  conformation,  with  its  channel  activity  simply  the 
expression  of  an  empty  state?  Three  protein  analogies  are  note¬ 
worthy:  prefoldin20,  protein  disulfide  isomerase  (PDI)21  and 
Salmonella  SicP22.  Neither  prefoldin  nor  PDI  require  ATP  for 
their  chaperone  activity.  The  prefoldin  hexamer  consists  of  six 
long  a-helical  coiled  coils  assembled  to  a  double  p-barrel  tem¬ 
plate.  The  distal  segments  of  the  coiled  coils  generate  the 
hydrophobic  surface  required  for  binding  of  unfolded  poly¬ 
peptides.  The  crystal  structure  of  prefoldin20  shows  similarity  to 
tire  long  helices  of  the  translocation  domain  of  BoNT  A9  and  B10 
{Z  score  =  5.0,  r.m.s.  devation  of  4.5  A  and  3.7  A,  respectively23). 
PDI,  whose  substrate-binding  cycle  is  regulated  by  its  redox  state, 
is  crucial  to  the  unfolding  of  cholera  toxin  and  for  its  transloca¬ 
tion  from  the  ER  lumen  into  the  cytosol  of  target  cells21. 
Numerous  bacterial  pathogens  deliver  proteins  (effectors)  into 
target  cells  with  the  assistance  of  specific  chaperones.  To  fit 
through  the  secretion  ‘channel’,  with  an  apparent  cut-off  diame¬ 
ter  of  30  A,  the  effector  is  kept  unfolded  by  interacting  with  the 
chaperone.  The  1.9  A  crystal  structure  of  the  chaperone-binding 
domain  of  the  Salmonella  effector  protein  SptP  in  complex  with 
its  cognate  chaperone  SicP  shows  the  unfolded  conformation  of 
SptP  threaded  around  three  chaperone  molecules22.  For 
BoNT  A  we  may  envision  a  chaperone  activity  driven  by  a  pH 
gradient  across  the  endosome  in  which  the  HC  channel  prevents 
tiie  aggregation  of  the  LC  in  the  acidic  vesicle  interior  via 
hydrophobic  interactions,  thereby  maintaining  the  unfolded 
conformation  during  translocation  and  releasing  the  LC  after  it 
refolds  at  the  neutral  cytosolic  pH.  In  this  process,  the  HC  chan¬ 
nel  is  ocduded  by  the  LC  during  transit  and  opens  after  comple¬ 
tion  of  translocation  and  release  of  cargo.  In  other  words,  the 
BoNT  A  HC-LC  complex  embedded  in  the  membrane  is  a 
dynamic  structural  device  that  prevents  aggregation  and  achieves 
translocation  of  the  LC.  An  analogous  mechanism  has  been 
invoked  for  the  translocation  of  the  catalytic  domain  of  diphthe¬ 
ria  toxin24.  There  is  a  dose  resemblance  of  the  ocduded  BoNT  A 
channel  with  the  translating  ribosome-protein-conducting 
channel  complex  (yeast  ER  sec61  complex):  the  cryo-EM  recon¬ 
struction  at  15.4  A  resolution  led  to  the  view  that  the  nascent 
polypeptide  chain  is  tightly  accommodated  within  the  channel, 
hindering  conductance25.  Given  a  constriction  of  -15  A,  the 
channel  may  allow  passage  of  only  a-helices25.  The  8  A  structure 
of  the  sec61 -related  bacterial  protein  translocation  complex 
SecYEG  reveals  a  deep  cavity  at  the  dimer  interface  that  is  occlud¬ 
ed  at  the  periplasmic  side  and  may  represent  the  dosed  transloca¬ 
tion  channel25.  Similarly,  unfolding  of  precursor  proteins  by  the 
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mitochondrial  import  system  is  required  for  translocation  across 
the  mitochondrial  Tim23  channel,  which  has  a  pore  diameter  of 
-13  A  (ref.  27).  For  chloroplasts,  the  peptide-translocating  chan¬ 
nels  of  the  outer  (TOC75)28  and  inner  (TicllO)29  membranes  dis¬ 
play  pore  diameters  of  -14  A  and  15  A,  respectively.  Protein 
conducting/translocating  channels  seem,  therefore,  to  constitute 
a  novel  mode  of  chaperone  action.  Our  findings  also  suggest  that 
the  BoNT  channel  may  represent  a  potential  target  for  interven¬ 
tion  to  attenuate  BoNT  neurotoxicity. 

Methods 

Reconstitution  of  BoNT  channels  in  lipid  bilayers.  Lipid  bilayers 
were  formed  at  25  ±  1  *C  by  apposition  of  two  monolayers30,  either  at 
the  tip  of  patch  pipets13  for  single  channel  recordings  (Figs.  1,  2)  or 
across  a  500  pm  aperture  on  a  teflon  partition  between  two  aqueous 
compartments  (250  pi)  for  LC  translocation  experiments12  (Figs.  3, 4). 
The  lipid  composition  was  optimized  to  enhance  the  binding  and 
insertion  of  BoNT  A  and  HC  and  consisted  of  (1,2-diphytanoyl- 
sn-glycero-3)-pho$phatidylethano!amine,  -phosphatidylcholine  and 
-phosphatidylserine  (Avanti  Polar  Lipids),  supplemented  with  GTIb 
gangliosides  (Fidia  Research  Laboratories)  at  a  weight  ratio 
2:2:1:0.125  for  the  four  lipids  used.  Acquisition  and  analysis  of  chan¬ 
nel  currents  were  essentially  as  described12*13.  Throughout  the  text, 
data  are  shown  as  mean  ±  s.d.  and  n  =  number  of  experiments. 

Native  BoNT  A  holotoxin,  LC  and  HC  Purified  native  BoNT  A, 
LC  A  and  HC  A  were  from  M.C  Goodnough  (University  of 
Wisconsin).  BoNT  A  nicking  was  estimated  by  SDS-PAGE  as  £99%. 

Recombinant  BoNT  A  LC  and  HC  BoNT  A  LC  clone  in  pBN3 
expression  vector  was  obtained  from  H.  Niemann.  The  LC  protein 
was  expressed  in  the  Escherichia  coli  strain  BL21  (Novagen)  and 
purified  essentially  as  described31.  The  recombinant  HC  clone  was 
obtained  from  B.R.  Singh  (University  of  Massachusetts)  and  sub¬ 
cloned  into  pCRT7/CT  TOPO  TA  vector  (Invitrogen).  The  HC  protein 
was  expressed  in  the  E.  coli  strain  BL21-CodonPlus(DE3)-RIL 
(Stratagene)  and  purified  from  the  6  M  guanidinium-CI-clarified 
inclusion  bodies  extract  by  nickel  chelate  chromatography  in  pres¬ 
ence  of  0.05%  dodecylphosphocholine  (Avanti  Polar  Lipids). 

Endopeptidase  activity  of  translocated  LC  In  bilayer  experi¬ 
ments.  Recombinant  SNAP-25  was  expressed  and  purified  as 
described32.  SNAP-25  truncation  was  evaluated  by  ELJSA.  N-terminal 


biotinylation  of  SNAP-25  was  accomplished  at  pH  72  using  the 
W-hydroxysuccinimide  ether  of  biotin  (biotin-XX-NHS,  Calbiochem). 
Cleavage  of  biotinylated  SNAP-25  (4  pM)  proceeded  for  24  h  at 
37  *C  ELISA  microtitration  plates  (Dynex  Technologies)  were  coated 
with  goat  polyclonal  antibodies  against  SNAP-25  C-terminal  pep¬ 
tide  (Santa  Cruz  Biotechnology).  Plates  were  blocked  with  biotin- 
free  1%  casein  colloid  buffer  (Research  Diagnostics).  Each  bilayer 
experimental  sample  was  tested  in  duplicate  and  in  parallel  to  con¬ 
trol  wells  on  the  same  plate,  after  incubation  for  1  h  at  37  °C. 
Streptavidin-polymerized  horseradish  peroxidase  (Poly-HRP-SAw) 
conjugate  (Research  Diagnostics)  was  applied  at  1:5,000  dilution. 
Absorbance  at  405  nm  was  determined  after  addition  of  ABTS  (2,2*- 
azino-bi$-3-ethylbenzothiazo!ine-6-sulfonic  add;  Boehringer 
Mannheim),  using  an  automatic  microplate  reader  (EL311;  Bio-Tek 
Instruments).  On  the  basis  of  the  standard  curve  (n  e  4),  the  ELISA 
detection  limit  of  <3x1 0-1*  M  was  established  as  the  BoNT  concen¬ 
tration  corresponding  to  3  s.d.  from  the  maximum  (100%) 
undeaved  SNAP-25  (Fig.  3d)53. 

CD  spectroscopy.  CD  measurements  were  taken  at  25  °C  on  an 
AVIV  202  CD  spectrometer  (Aviv  Instruments)  using  a  0.1  cm  path 
length  cell  and  analyzed  as  described34.  For  all  spectra,  the  average 
of  three  scans  was  obtained.  To  correlate  the  CD  results  with  the 
bilayers  measurements,  we  used  100  mM  ICO,  50  mM  potassium 
phosphate  buffered  to  pH  7.0,  5.0  or  4.5.  These  buffers  were  sup¬ 
plemented  with  1  mM  p-mercaptoethanol  and  0.05%  dodecylphos¬ 
phocholine  to  prevent  aggregation,  which  increases  with 
acidification.  The  differences  between  our  CD  results  and  other 
reports35  likely  arise  from  the  different  experimental  conditions 
used.  Li  and  Singh35  used  50  mM  Nad  buffered  with  50  mM  Tris  for 
the  pH  7.0  measurements  and  50  mM  NaCI  buffered  with  20  mM  Tris 
for  the  pH  4.7  measurements,  and  their  measurements  were  con¬ 
ducted  at  4  °C 
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ABSTRACT  Complex  conformational  changes  influence  and  regulate  the  dynamics  of  ion  channels.  Such  conformational 
changes  are  stochastic  and  often  inhomogeneous,  which  makes  it  extremely  difficult,  if  not  impossible,  to  characterize  them  by 
ensemble-averaged  experiments  or  by  single-channel  recordings  of  the  electric  current  that  report  the  open-closed  events  but 
do  not  specifically  probe  the  associated  conformational  changes.  Here,  we  report  our  studies  on  ion  channel  conformational 
changes  using  a  new  approach,  patch-clamp  fluorescence  microscopy,  which  simultaneously  combines  single-molecule 
fluorescence  spectroscopy  and  single-channel  current  recordings  to  probe  the  open-closed  transitions  and  the  conformational 
dynamics  of  individual  ion  channels.  We  demonstrate  patch-clamp  fluorescence  microscopy  by  measuring  gramicidin  ion 
channel  conformational  changes  in  a  lipid  bilayer  formed  at  a  patch-clamp  micropipette  tip  under  a  buffer  solution.  By  measuring 
single-pair  fluorescence  resonance  energy  transfer  and  fluorescence  self-quenching  from  dye-labeled  gramicidin  channels,  we 
observed  that  the  efficiency  of  single-pair  fluorescence  resonance  energy  transfer  and  self-quenching  is  widely  distributed, 
which  reflects  a  broad  distribution  of  conformations.  Our  results  strongly  suggest  a  hitherto  undetectable  correlation  between 
the  multiple  conformational  states  of  the  gramicidin  channel  and  its  closed  and  open  states  in  a  lipid  bilayer. 


Subtle  conformational  changes  of  ion  channels  play  an 
important  role  in  regulating  the  function  and  dynamics  of 
ion  channels  (Karlin,  2002;  Madden,  2002;  Choe,  2002). 
Channel  conductance  recording  analysis  and  kinetic  model¬ 
ing  (Cross  et  al.,  1999;  Andersen  et  al.,  1999),  together  with 
site-directed  mutagenesis,  protein  modifications,  and  static 
structural  analysis  provide  extensive  insights  into  the 
conformations  and  the  conformational  rearrangements  as¬ 
sociated  with  channel  activity  (Karlin,  2002;  Madden,  2002; 
Choe,  2002;  Glauner  et  al.,  1999;  Cha  et  al.,  1999). 
However,  the  conformational  dynamics  underlying  the 
mechanisms  of  ion  channel  action  are  still  largely  unknown 
due  to  lack  of  direct  measurements.  Mechanistic  character¬ 
ization  of  the  conformational  changes  has  been  primarily 
based  on  kinetic  model  analyses  of  the  electric  current 
trajectories  obtained  by  single-channel  and  whole-cell  patch- 
clamp  recordings  that  are  not  sensitive  to  subtle  conforma¬ 
tional  intermediate  states,  including  electrically  undetectable 
“silent”  conformational  states.  Protein  dynamics  may  in¬ 
volve  significant  static  and  dynamic  inhomogeneities  (Lu 
et  al.,  1998;  Xie  and  Trautman,  1998).  For  instance,  ion 
channel  dynamics  can  be  modulated  by  inhomogeneous 
local  membranes  and  a  fluctuating  environment.  These  static 
and  dynamic  inhomogeneous  conformational  dynamics  are 
extremely  difficult  to  study  by  conventional  experimental 
approaches. 
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Recently,  our  group  (Ott  et  al.,  2001,  2002;  Harms  et  al., 
2002,  2003)  and  others  have  made  significant  progress  in 
developing  new  approaches  to  studying  ion  channels  using 
both  fluorescence  imaging  and  electric  patch  recording. 
Utilizing  a  micro-pinhole  patch  technique,  Yanagida  and 
coworkers  (Ide  and  Yanagida,  1999;  Ide  et  al.,  2002)  first 
reported  fluorescence  images  acquired  before  or  after 
a  single-channel  current  measurement  on  the  same  sample 
lipid  bilayer.  Woolley,  Schuetz,  and  their  co-workers 
(Lougheed  et  al.,  2001;  Borisenko  et  al.,  2003)  have  ele¬ 
gantly  demonstrated  the  feasibility  of  combining  single¬ 
molecule  fluorescence  imaging  and  electric  recording  on 
gramicidin  ion  channels.  Isacoff  and  co-workers  (Sonnleit- 
ner  et  al.,  2002)  have  reported  a  single-molecule  imaging  and 
whole-cell  patch  recording  study  of  structural  rearrangement 
of  voltage-gated  Shaker  K+  channels  in  the  plasma 
membrane  of  living  cells.  A  related  experimental  effort  by 
Schmidt  and  his  co-workers  (Harms  et  al.,  2001)  has  dem¬ 
onstrated  wide-field  fluorescence  microscopy  imaging  of  the 
diffusion  and  aggregation  of  yellow-fluorescent  protein- 
labeled  L-type  Ca2+  channels  in  living  cells.  Another  related 
experiment  by  Zagotta  and  his  co-workers  (Zheng  and 
Zagotta,  2000)  has  studied  gating  rearrangements  of  cyclic 
nucleotide-gated  channel  proteins  by  combined  fluorescence 
imaging  and  patch-clamp  recording  at  multiple-channel- 
averaged  measurements.  These  research  efforts  have  begun 
to  provide  a  new  paradigm  for  studying  ion  channel  confor¬ 
mational  dynamics  and  mechanisms  and  for  obtaining  in¬ 
formation  not  obtainable  by  the  conventional  patch  current 
recording  measurements,  a  primary  experimental  approach 
for  decades. 

Here,  we  demonstrate  a  new  approach,  patch-clamp 
fluorescence  microscopy  (PCFM),  by  combining  single¬ 
molecule  fluorescence  imaging  of  fluorescence  resonant 
energy  transfer  (FRET)  and  fluorescence  self-quenching 


Single  Ion  Channel  Dynamics 

with  a  single-channel  electric  current  patch  recording.  Using 
PCFM,  we  were  able  to  correlate  in  real-time  the  single¬ 
channel  open-closed  kinetics  with  simultaneous  changes  in 
optical  signals  associated  with  conformational  changes  in 
a  gramicidin  ion  channel  in  a  lipid  bilayer.  By  utilizing 
a  typical  patch-clamp  technique  in  this  new  approach,  the 
application  of  PCFM  for  studying  ion  channel  dynamics  in 
living  cells  (On*  et  al.,  2002)  is  now  feasible.  A  two-state 
kinetic  model  has  been  widely  referenced  to  illustrate  the 
gramicidin  ion  channel  kinetic  behavior,  i.e.,  the  association 
and  dissociation  of  the  gramicidin  dimer  corresponding  to 
the  channel’s  open  and  closed  states.  However,  the  nature  of 
the  conformational  changes  and  the  dynamics  that  regulate 
channel  activity  are  still  largely  unclear,  since  there  are  no 
methods  of  providing  a  direct  real-time  single-molecule 
measurement  of  both  ion  channel  electric  current  and 
conformational  changes.  Our  new  experimental  results,  not 
obtainable  by  using  ensemble-averaged  methods  or  single¬ 
channel  patch  electric  recording  alone,  suggest  the  oc¬ 
currence  of  multiple  intermediate  conformational  states 
underlying  the  gramicidin  channel  dynamics. 

MATERIALS  AND  METHODS 

Linear  gramicidin  produced  naturally  from  the  bacterium  Bacillus  brevis 
(Hladky  and  Haydon,  1970)  is  a  typical  “model”  ion  channel  (Koeppe  and 
Andersen,  1996;  Cross  et  al.,  1999;  Andersen  et  al.,  1999),  and  its 
conductive  open  state  is  presumably  a  dimer  (Veatch  et  al.,  1975;  Koeppe 
and  Andersen,  1996;  Cross  et  al.,  1999;  Andersen  et  al.,  1999)  of  two  15- 
amino  acid  peptides  forming  /3-helices  that  meet  head-to-head  at  the 
N-termini  in  the  interior  of  a  lipid  bilayer  (Koeppe  and  Andersen,  1996; 
Cross  et  al.,  1999;  Andersen  et  al.,  1999;  Roux  and  Kaiplus,  1994).  Because 
of  the  available  extensive  knowledge  on  the  gramicidin  ion  channels  in 
lipid  bilayers,  we  chose  this  system  to  demonstrate  our  combined  single¬ 
molecule  approach  of  PCFM. 
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Fig.  1  shows  the  experimental  setup  for  the  PCFM  single-channel  electric 
current  recording  and  fluorescence  imaging.  We  recorded  the  single-channel 
current  of  dye-labeled  gramicidin  dimers  that  were  incorporated  in  lipid 
bilayers  formed  at  the  tip  of  a  patch  pipette  (Montal  and  Mueller,  1972;  see 
Fig.  1,  inset)  as  we  simultaneously  imaged  the  pipette  tip  using  wide-field 
single-molecule  fluorescence  microscopy  (Schmidt  et  al.,  1996;  Ma  et  al., 
2000;  Han  et  al.,  2001;  Bartko  et  al.,  2002;  Seisenberger  et  al.,  2001; 
Deschenes  and  Vanden  Bout,  2001;  Cognet  et  al.,  2000). 


Single-molecule  fluorescence  microscope 

A  microscope  (Nikon  Diaphot  300,  Nikon,  Kanagawa,  Japan)  equipped 
with  a  60 X  objective  (N.A.  =  1.4)  and  4X  relay  lens  was  used  in  this 
work.  In  all  experiments,  the  illumination  intensity  was  set  to  5  kW/cm2 
for  tetramethylrhodamine  (TMR)-gramicidin  at  514  nm  (Ar-ion  laser. 
Molecular  Probes,  Eugene,  OR;  MWK  Laser,  Riverside,  CA)  or  532  nm 
(Verdi  10  W  diode  laser,  Coherent,  Palo  Alto,  CA),  and  1.5  kW/cm2  for 
Cy5 -gramicidin  at  632  nm  (HeNe  laser,  Uniphase  or  dye  laser,  CR-5999, 
Coherent).  The  polarization  of  the  excitation  light  was  controlled  by  using 
a  A/4  waveplate  (Meadowlark  Optics,  Frederick,  CO).  Filter  combinations 
for  TMR  consisted  of  DCLP555,  HQ625/100  (Chroma  Technology  Corp., 
Brattleboro,  VT)  and  OG550  (Schott,  Duiyea,  PA),  and  for  Cy-5  of 
DCLP645,  HQ690/90,  and  RG645  (Schott).  The  single-pair  fluorescence 
resonant  energy  transfer  (spFRET)  measurements  were  done  using 
a  custom-built  dual-dichroic  beam-splitter  (a  combination  of  a  DCLP555 
and  DCLP645,  Chroma  Technology  Corp.)  and  band-pass  filter  in 
addition  to  either  OG550  or  RG645,  permitting  the  detection  of  single 
fluorophores  by  a  nitrogen-cooled  CCD  camera  (Spec- 10  1340  X  400B  or 
700B,  Roper  Scientific,  Trenton,  NJ).  The  microscope  system  total 
detection  efficiency  was  10  ±  2%  for  TMR  and  8  ±  2%  for  Cy5.  We 
estimated  the  total  detection  efficiency  from  the  collection  efficiency  of 
the  objective,  the  spectral  throughput  of  the  filter  set,  and  the  spectral 
detection  efficiency  of  the  CCD  camera.  We  also  subtracted  the  estimated 
reflections  for  each  glass  surface  in  the  optical  beam  path.  This  device 
ensured  the  toggling  of  the  excitation  wavelength  between  532  nm  and 
632  nm,  and  the  detection  of  the  emission  at  8  ±  2%  efficiency  for  TMR 
and  5  ±  2%  efficiency  for  Cy5  with  the  FRET  filters.  The  images  were 
taken  using  5-ms  exposures  at  1-10  Hz. 
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FIGURE  1  Experimental  apparatus  and  mea¬ 
surement  configuration  of  PCFM.  Lipid  bilayers 
were  formed  at  the  tip  of  a  patch  pipette  by 
apposition  of  two  monolayers  of  mixed  lipids  (4:1 
dPhPE:DPhPC).  Dye-labeled  gramicidin  mono¬ 
mers  are  introduced  from  both  sides  of  the  bilayer, 
and  single-channel  electric  currents  are  recorded  by 
using  typical  patch-clamp  instrumentation.  The 
laser  excitation  was  focused  through  a  microscope 
objective  on  the  pipette  tip,  and  the  fluorescence 
was  imaged  on  a  CCD  camera. 
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Lipid  bilayer  formation  at  the  tip  of  a 
patch-pipette 

Lipid  bilayers  were  formed  at  the  tip  of  patch  pipettes  by  the  apposition  of 
two  monolayers  (Montal  and  Mueller,  1972),  using  4:1  diphytanoylphos- 
phatidylethanolamine  (dPhPE):diphytanoylphosphatidylcholine  (dPhPC) 
(Avanti,  Alabaster,  AL).  Gramicidin  dissolved  in  ethanol  (1CT9  M)  was 
added  to  the  aqueous  subphase  (1  M  KCI,  1  mM  CaCl2>  and  5  mM  Hepes  of 
pH  7.5)  after  bilayer  formation  (Fig.  1,  inset).  Gramicidin  monomers  were 
labeled  with  either  TMR  or  Cy5.  Homodimers  of  TMR-labeled  gramicidin 
were  assembled  by  adding  gramicidin-TMR  monomers  to  both  sides  of  the 
membrane;  gramicidin-TMR-gramicidin*Cy5  heterodimers  were  generated 
by  introducing  TMR-gramicidin  and  Cy5-gramicidin  from  opposite  sides  of 
the  bilayer  (O’Connell  et  al.,  1990).  The  latency  period  to  form  an  active 
channel  is  typically  a  few  minutes  after  gramicidin  insertion  into  the  bilayer 
and  presumably  arises  from  the  slow  diffusion  of  individual  gramicidin 
monomers  to  form  a  dimer. 


Single-channel  electric  current  patch  recording 

Single-channel  electric  currents  elicited  at  constant  applied  voltage  (—100- 
+  100  mV)  were  recorded  using  an  EPC-9  amplifier  (HEKA  Electronics, 
Lambrecht,  Germany).  The  recorded  channel  electric  signal  was  sampled  at 
4  kHz  and  filtered  at  3  kHz.  Control  experiments  using  unlabeled  gramicidin 
A  revealed  distributions  of  current  amplitude  and  open/closed  dwell-time 
similar  to  those  of  the  dye-labeled  gramicidin.  The  dye-labeling  of  the 
gramicidin  C-terminal  did  not  significantly  alter  the  gramicidin  ion  channel 
activities  (Lougheed  et  al.,  2001;  Borisenko  et  al.,  2003).  The  presence  of 
a  single  channel  at  the  patch  was  confirmed  by  recorded  channel-current 
trajectories  showing  typical  single-channel  open-closed  behavior  (O’Con¬ 
nell  et  al.,  1990;  Sawyer  et  al.,  1989;  Veatch  et  al.,  1975).  The  following 
criteria  were  applied  to  ensure  that  the  recorded  channel  activity  was  from 
a  single  ion  channel  and  not  due  to  contamination  or  membrane  instability: 
1),  no  channel  activity  was  detected  before  adding  gramicidin;  2),  only 
bursts  of  channel  activity  were  considered;  3),  a  high  membrane  resistance 
(at  the  giga-ohm  range)  leading  to  a  high  signal-to-noise  ratio  (S/N)  was 
achieved;  and  4),  membrane  breakdown  occurred  when  300-500  mV 
polarization  voltage  was  applied  to  the  patch.  These  criteria  are  also 
practiced  commonly  in  other  labs  (O’Connell  et  al.,  1990;  Sawyer  et  al., 
1989;  Veatch  et  al.,  1975;  Sigworth  et  al.,  1987)  to  ensure  a  single 
gramicidin  channel  patch  recording.  About  10%  of  all  the  attempts  were 
successful  in  incorporating  a  single  gramicidin  channel  at  the  patch; 
membranes  without  channel  activity  or  with  more  than  one  channel  were  not 
pursued. 


Dye  labeling  of  gramicidin  monomers 

Gramicidin  C  with  Gly  and  Lys  residues  attached  to  the  C-terminal 
(synthesized  by  Genemed  Synthesis,  Inc.)  was  used  in  this  work.  TMR  or 
Cy5  HNS  esters  were  covalently  bonded  by  acylation  to  the  e-amino  group 
of  the  Lys  residue.  The  full  sequence  of  the  gramicidin  used  for  the 
fluorescence  labeling  was  Formyl  -  V  -  G  -  A  -  (d)L  -  A  -  (d)V  -  V  -  (d)V  -  W  - 
(d)L  -  Y  -  (d)L  -  W  -  (d)L  -  W  -  G  -  K-NHCH20H.  Multiple  HPLC 
purifications  (C8  column  with  80%  methanol  solvent)  ensured  the  purity  of 
the  dye-labeled  gramicidin  samples,  which  was  further  confirmed  by 
MALDI-TOF  mass  spectroscopy  (molecular  weight  of  gramicidin  C-TMR, 
calculated  =  2412  and  measured  =  2415;  molecular  weight  of  gramicidin 
C-Cy5,  calculated  =  2637  and  measured  =  2641). 


Fluorescence  imaging  data  analyses 

The  fluorescence  background  was  subtracted  by  averaging  photon  counts  for 
each  pixel  or  by  fitting  the  background  to  a  wide  2-D  Gaussian.  The 


individual  fluorescent  “hot  spots”  were  identified  by  the  maxima  of  a  2-D 
Gaussian-correlation  filter  (Schmidt  et  al.,  1996).  The  region  around  these 
maxima  was  fitted  to  a  2-D  Gaussian,  reporting  the  intensity,  width,  position, 
and  errors  of  the  least  square  fit  (Schmidt  et  al.,  1995,  1996;  Press  et  al., 
1990). 

Single-molecule  fluorescence 
control  experiments 

Single-molecule  imaging  was  conducted  of  dye-labeled  gramicidin  in 
Langmuir-Blodgett  (LB)  bilayers  and  spin-coated  over  polymethyl  methya- 
crylate  (PMMA)  polymer-coated  cover  slips.  The  single-molecule  fluores¬ 
cence  imaging  signal  was  confirmed  by  several  typical  verification 
measurements,  including  observations  of  single-step  photobleaching  events 
(Fig.  2  ay  inset)  and  of  linearly  polarized  dipoles  (data  not  shown).  The 
dPhPE: dPhPC  4:1  was  doped  with  1(T8  mol  TMR-gramicidin  C  or  Cy5- 
gramicidin  C  per  mol  lipid  and  were  deposited  at  22°C  and  28  dynes/cm  on 
glass  slides  using  LB  deposition.  Single-molecule  fluorescence  imaging 
measurements,  using  essentially  the  same  experimental  parameters  and 
configurations  that  were  used  to  image  the  bilayers  containing  individual 
dye-labeled  gramicidin  molecules  at  the  tip,  provided  a  quantitative 
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FIGURE  2  (a)  Histogram  of  single  hot-spot  fluorescence  intensities 
deduced  from  2-D  Gaussian  fitting  to  each  spot  imaged  at  the  patch-clamp 
tip.  Inset:  intensity  trajectory  from  a  “hot  spot”  image  of  a  single  TMR- 
gramicidin  molecule,  (b)  The  probability  density  function  ( normalized ) 
deduced  from  the  histogram  in  a ,  and  from  the  histogram  of  hot-spot 
intensities  on  glass-supported  lipid  bilayers  (< dotted  line)  as  a  control.  Our 
control  experiments  and  those  of  others  (Lougheed  et  al.,  2001)  show  no 
observable  perturbation  by  dye  labeling  on  the  gramicidin  channel  activity. 
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reference  for  the  single-molecule  photon  count  rates  under  the  wide-field 
imaging  configuration  (Fig.  2  a)  in  our  experiments,  which  is  a  typical 
calibration  commonly  practiced  in  single-molecule  imaging  experiments 
(Schmidt  et  al.,  1995, 1996;  Seisenberger  et  al.,  2001;  Harms  et  al.,  2001). 
The  single-molecule  images  were  analyzed  by  the  least-square  2-D  Gaussian 
fitting,  as  discussed  above. 

The  probability  density  function  (PDF)  of  the  photon  counts, 


PDF(x)  =  £ 


A/iV2^ 


2A/f  / 


is  constructed  using  parameters  from  the  single-molecule  image  photon 
counts,  where  I„  A /*,  and  represent  the  photon  counts,  its  standard 
deviation,  and  its  relative  weight.  The  PDF  statistically  evaluates  the 
imaging  photon  counts  with  their  standard  deviations  and  occurrences  and 
provides  a  more  reliable  distribution  than  an  occurrence  histogram  does 
(Schmidt  et  al.,  1995).  Under  our  set-up  and  experimental  conditions,  die 
average  photon  count  detection  rate  for  the  TMR-gramicidin  C  monomers 
was  68  ±  5  counts/ms  with  a  PDF  peak  of  45  ±  5  counts/ms  and  PDF  width 
of  32  counts/ms.  For  Cy5-gramicidin  C,  the  average  photon  count  was  46  ± 
5  counts/ms  with  a  PDF  peak  of  36  ±  5  counts/ms  and  width  of  30  counts/ 
ms.  The  width  of  the  PDF  distributions  indicates  that  the  level  of  signal 
fluctuations  is  due  primarily  to  shot  noise  and  to  the  fluorescence  quantum 
yield  distribution  of  the  individual  dye  molecules  in  a  spatially  inhomo¬ 
geneous  local  environment  (Lu  and  Xie,  1997;  Xie  and  Trautman,  1998). 

There  are  inevitable  differences  in  the  fluorescence  photon  detection 
efficiency  and  the  signal-to-noise  ratio  of  imaging  single  molecules  on  the 
LB  layer  or  the  PMMA  polymer,  as  in  a  control  experiment,  and  in  imag¬ 
ing  single  molecules  at  a  patch-clamp  pipette  tip,  as  in  a  real  PCFM 
measurement  (Fig.  2  b).  There  are  at  least  three  reasons  for  these  differences: 
1),  There  are  different  configurations  for  the  pipette  tip  under  a  buffer 
solution  in  the  PCFM  measurement  and  the  cover  slip  with  individual  dye- 
labeled  gramicidin  molecules  in  LB  layers  under  water  or  spin-coated 
individual  dye-labeled  gramicidin  molecules  on  PMMA  polymer  films. 
These  intrinsic  differences  in  the  experimental  configuration  apparently 
cause  different  signals  and  noise  levels  associated  with  a  change  in  the 
refractive  index  and  light-scattering  conditions.  2),  There  are  different 
fluorescence  quantum  efficiencies  under  different  conditions  for  the  LB 
bilayer,  a  polymer,  and  the  lipid  bilayer  at  a  pipette  tip  under  a  buffer.  3), 
There  are  different  focal  planes  of  the  wide-field  laser  excitation  and 
fluorescence  collection,  which  are  at  the  upper  surface  of  the  cover  slip  for 
the  control  experiment  but  at  —70  jjl m  above  the  cover  slip  in  the  buffer  for 
imaging  at  the  pipette  tip. 

The  single-molecule  fluorescence  signal  was  distinguished  from  the 
background  at  the  pipette  tip  based  on  characteristic  single-molecule 
spectroscopic  signatures:  fluorescence  blinking,  single-step  photobleach- 
ing,  specific  polarization  orientation,  and,  most  importantly,  simultaneous 
single-channel  ion  current  recording.  The  last  simultaneous  verification  was 
that  only  a  single-channel  ion  current  associated  with  a  single  ion  channel 
in  the  patch  was  observed  during  fluorescence  imaging.  The  signal-to-noise 
ratio  was  typically  within  a  range  of  5-10  in  the  control  LB  bilayer 
imaging  at  the  LB  bilayer  or  PMMA  polymer  film,  and  3-6  in  the  PCFM 
imaging  at  the  pipette  tips. 


RESULTS  AND  DISCUSSION 

Single-molecule  spectroscopy  is  powerful  in  resolving 
complex  and  inhomogeneous  dynamics  and  kinetics  of 
protein  systems.  However,  the  fluctuation  of  most  of  the 
measured  parameters  is  intrinsically  pertinent  to  single¬ 
molecule  measurements  (Zwanzig,  1990;  Lu  and  Xie,  1997; 
Xie  and  Trautman,  1998;  Moemer,  2002).  For  example, 
single-molecule  fluorescence  has  spectral  and  intensity 
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fluctuations,  FRET  efficiency  fluctuations,  and  lifetime 
fluctuations.  These  fluctuations  are  not  observable  in 
a  typical  ensemble-averaged  experiment,  and  in  most  cases 
they  are  associated  with  single-molecule  conformations  and 
local  environments.  A  reliable  conclusion  should  typically 
be  drawn  from  sound  statistical  analyses  of  a  sufficient 
number  of  measurements  of  single-molecule  events. 
Although  the  fluctuations  give  a  larger  standard  deviation 
in  the  data  of  a  single-molecule  measurement  compared  to 
that  of  an  ensemble-averaged  measurement,  the  information 
obtained  from  single-molecule  experiments  has  not  been 
obtainable  from  ensemble-averaged  experiments. 

Heterodimers  of  gramicidin  statistically  correlate 
with  maximum  spFRET  at  an  open  state 

Using  PCFM,  we  conducted  colocalization  and  spFRET 
imaging  experiments  simultaneously  with  single-channel 
patch  recording.  Asymmetric  incorporation  of  Cy 5 -gramici¬ 
din  C  and  TMR-gramicidin  C  led  to  the  formation  of  single 
channels  tagged  with  two  different  dyes,  the  donor  and 
acceptor  of  the  fluorescence  energy  transfer.  Colocalization 
and  spFRET  within  a  single  gramicidin  heterodimer  were 
probed  by  dual-color  excitation  and  donor-acceptor  two- 
channel  fluorescence  imaging,  using  an  instrument  config¬ 
uration  similar  to  one  reported  previously  (Cognet  et  al., 
2000).  Fig.  3  a  shows  two  fluorescence  images  taken  at  the 
tip  with  514-nm  (image  1)  and  632-nm  (image  2)  excitation 
wavelengths,  using  a  filter  that  allows  imaging  of  both  Cy5 
and  TMR  emissions,  both  of  which  are  clearly  discemable  in 
Fig.  3  a.  From  the  2-D  Gaussian  fitting  of  the  colocalized 
signals,  the  two  dye  molecules  were  found  to  be  within  the 
fluorescence  imaging  diffraction-limited  spot  (~300-nm 
width),  which  suggests  that  a  pair  of  colocalized  TMR- 
gramicidin  and  Cy5-gramicidin  monomers  was  present  at 
the  tip  when  the  active  single-channel  current  signal  was 
observed  (Fig.  3  a ).  The  colocalization  imaging  measure¬ 
ment  is  important  and  was  used  periodically  to  check  the 
coexistence  of  the  donor  and  acceptor  molecules  and  to 
differentiate  the  situations  of  non-FRET  from  acceptor 
photobleaching  events.  To  further  characterize  the  hetero¬ 
dimerization  states  of  TMR-gramicidin  and  Cy5-gramicidin, 
spFRET  images  were  collected  using  a  645 -nm  long-pass 
filter,  allowing  only  the  fluorescence  of  acceptor  Cy5  to  pass. 
Toggling  between  514-nm  and  632-nm  excitation  wave¬ 
lengths  enabled  the  measurement  of  spFRET  and  verification 
of  the  coexistence  of  the  single  acceptor  (Fig.  3  b).  The 
spFRET  that  occurred  during  the  open  states  of  the  channel  is 
shown  in  Fig.  3  b.  By  counting  images  correlated  with  the 
electric  signal  that  reflected  the  open-closed  activity  of  single 
channels  at  the  patch,  we  found  that  —70%  of  channel¬ 
opening  events  were  correlated  with  spFRET.  Estimating 
that  ~25%  of  the  molecules  undergo  blinking  events,  we 
attributed  the  small  fraction  of  non-FRET  to  acceptor 
photobleaching.  Interestingly,  we  also  observed  that  spFRET 
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FIGURE  3  PCFM  measurement  of  single  TMR/Cy5- 
gramicidin  heterodimer  ion  channels.  Green  and  red  arrows 
marie  the  timing  of  the  toggled  excitation  at  514  nm  and 
632  nm,  respectively,  (a)  Colocalization  of  two  gramicidin 
monomers  labeled  with  two  different  dyes,  TMR  and  Cy5. 
Single-channel  current  trajectory  recorded  at  100-mV 
polarization  potential,  with  arrows  indicating  the  point  of 
exposure  of  two  consecutively  recorded  frames.  The  first 
frame  was  taken  at  514-nm  excitation  (with  the  OG550 
long-pass  filter),  indicating  single-molecule  fluorescence 
of  TMR.  The  second  frame  was  taken  at  632-nm 
excitation,  indicating  single-molecule  fluorescence  of 
Cy5.  From  a  2-D  Gaussian  fitting,  the  molecules  were 
found  to  be  within  30  nm  of  each  other  (the  minimal 
positional  accuracy  allowed  by  the  measurement  system). 
“C”  and  “0”  denote  the  current  corresponding  to  closed 
and  open  states,  respectively.  ( b )  spFRET  within  a  pair  of 
gramicidin  monomers  labeled  with  two  different  dyes, 
TMR  and  Cy5.  Single-channel  current  trajectory  (upper 
panel)  with  four  arrows  indicating  the  exposure  of  four 
consecutively  acquired  frames  (lower panel).  Frames  1  and 
3  were  taken  at  514-nm  excitation  with  a  red-emission 
long-pass  filter  (RG645).  Frames  2  and  4  were  taken  at 
632-nm  excitation  with  a  red-emission  long-pass  filter 
(RG645).  Frame  1  indicates  that  FRET  occurred  between 
the  TMR/Cy5  spFRET  pair  when  the  gramicidin  channel 
was  at  an  open  state.  Frame  2  verifies  that  a  single  Cy5 
molecule  was  present.  Frame  3  indicates  that  no  spFRET 
occurred  when  the  channel  was  at  a  closed  state.  Frame  4 
verifies  that  the  Cy5  was  still  present.  Subsequent  images 
taken  without  the  red  long-pass  filter  (RG645)  did  show 
colocalization  as  shown  in  a.  (c)  FRET  efficiency 
distribution  of  “channel  open”  TMR/Cy5  heterodimers 
( shaded  bars)  with  mean  of  0.59,  width  of  0.21,  and 
“channel  closed”  (open  bars)  with  a  maximum  value  of 
0.57.  An  asterisk  is  inserted  to  indicate  the  70%  of  channels 
showing  no  appreciable  FRET  (below  0.3)  in  the  channel- 
closed  state  (or  30%  of  channels  in  the  channel-open  state). 


FRET  efficiency 


e 


occurred  in  >—10%  of  the  cases  where  the  channels  were 
closed.  We  have  further  verified  the  FRET  anticorrelation  of 
the  donor-acceptor  fluorescence  intensities  by  photon  time- 
stamping,  time-correlated  single-photon  counting  (Hu  and 
Lu,  2003)  and  proved  the  presence  of  FRET  (Harms  et  al., 
2003;  G.  S.  Harms,  G.  Orr  and  H.  P.  Lu,  unpublished 
results). 

A  broad  distribution  of  FRET  efficiency  (Ha  et  al.,  1999; 
Brasselet  et  al.,  2000;  Cognet  et  al.,  2000), 

F—(  1  +  ^CyS  frey 5  ^TMr\ 

V  %MR  ^TMR  IcyS  ) 

was  found  at  gramicidin  heterodimers,  and  the  efficiency  was 
evaluated  by  the  activity  of  the  channels  at  the  time  that  the 
images  were  taken  (Fig.  3  c).  7™ R  and  /Cy5  are  the  total 
photon  counts  for  single-molecule  images  of  TMR  and  Cy5, 
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7]  is  detection  efficiency,  and  <3>  is  the  quantum  yield.  We 
used  the  values  of  7]tmr  and  Vcys  discussed  above  as  well  as 
^tmr  =  ^Cy5  =  0.28  (Cognet  et  al.,  2000).  When  the 
channels  were  open,  the  Gaussian-like  distribution  of  the 
efficiency  showed  a  mean  value  of  0.59  ±  0.05  with  a  full- 
width  half-maximum  (FWHM)  of  0.21.  The  standard 
deviation  of  the  F  was  estimated  to  be  ±0.05  based  on 
shot-noise  error  propagation  analysis  (Press  etal.,  1990).  Our 
control  experiment  of  polarization-dependent  imaging 
(Harms  et  al.,  2003)  indicated  that  the  rotational  diffusion 
of  the  TMR  and  Cy5  dye  molecules  had  a  much  faster  rate 
than  the  imaging  rate,  so  that  the  orientation  effect  of  the  dye 
dipoles  on  the  FRET  efficiency  was  averaged  out  during  the 
5-ms  imaging  collection  time  for  each  frame.  The  mean 
spFRET  efficiency  and  the  FWHM  of  0.21  correspond  to 
a  mean  donor-acceptor  distance  (d)  of  —56  A  and  a  range  of 
—52  A  <  d  <  —63  A,  respectively.  However,  —30%  of 
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the  images  obtained  from  the  channel  open  states  showed  no 
detectable  FRET  efficiency,  mostly  because  of  photobleach- 
ing  and  spectral  fluctuation.  We  cannot  differentiate  the  data 
between  no  FRET  and  low-efficiency  FRET  below  0.3,  due 
to  the  experimental  limitation  of  the  signal-to-noise  ratio 
of  the  single-molecule  PCFM  imaging.  The  distribution  of 
FRET  efficiency  of  the  closed  state  of  single  channels  had 
a  mean  of  0.42  ±  0.05  (excluding  the  low  signal  data)  and 
extended  to  a  maximum  efficiency  of  0.57.  Therefore,  the 
difference  in  the  FRET  between  open  and  closed  states  was 
statistically  significant  and  beyond  the  measurement  noise. 
The  broad  FRET  distributions  of  both  open  and  closed  states 
of  the  single  gramicidin  channels  suggest  that  the  channels 
could  be  in  closed  states  even  when  the  two  gramicidin 
monomers  are  still  in  an  intermediate  dimerized  state  (Fig.  4, 
a  and  b).  We  postulate  that  the  partially  dissociated  dimer 
may  have  had  a  distorted  ion  channel  pore,  creating 
a  significant  change  in  the  channel  electric  current,  though 
the  physical  distance  between  the  two  monomers  did  not 
significantly  change  (Fig.  4  b).  It  has  been  suggested  that 


FIGURE  4  Molecular  stnictures  (top  and  side  views)  of  the  possible  open 
{a)  and  closed  (b)  conformational  states  of  a  gramicidin  dimer  in  a  lipid 
bilayer.  The  molecular  structures  were  calculated  by  molecular  dynamics 
simulation  using  Discovery  and  Insight  O  software.  It  is  evident  that  slight 
rotational  distortion  along  the  dihedral  angle  can  significantly  narrow  the 
channel  pore  size  and  separate  the  dye  probe  pair. 


there  are  six  hydrogen  bonds  at  the  N-terminal  head-to-head 
cross-linking  the  dimer  (Szabo  and  Urry,  1978;  Koeppe  and 
Andersen,  1996).  It  is  possible  that  one  or  more  broken 
hydrogen  bonds  can  give  different  conformations  of  the  pore 
structure  of  a  gramicidin  dimer.  In  Fig.  4,  we  show  that 
the  energy-minimized  structure  of  gramicidin  heterodimer 
FRET  pair  (represented  by  using  TMR  molecular  structures) 
has  the  minimum  donor-acceptor  distance  of  ~50-52  A.  The 
relative  dihedral  angle  between  the  two  monomers  at 
different  conformations  of  the  conducting  dimer  (Fig.  4) 
can  give  a  distance  change  up  to  6  A,  which  can  result  in 
a  0.15  change  in  FRET  efficiency.  We  note  that  the 
calculated  spFRET  efficiency  is  only  reliable  in  the  range 
of  0.3  <  F  <  0.7  due  to  the  measurement  noise  associated 
with  the  effects  of  shot  noise  and  photophysical  dark  states. 
Nevertheless,  spFRET  imaging  results  (Fig.  3  c)  provide 
unambiguous  evidence  that  the  gramicidin  channel  dynamics 
involve  multiple  open  and  closed  states  associated  with 
different  conformations  of  the  gramicidin  dimers. 

The  spFRET  imaging  measurement  using  toggled  dual¬ 
color  excitation  and  simultaneous  correlation  with  the  single 
gramicidin  ion  channel  patch  electric  recording  is  a  novel 
approach,  and  the  results  of  the  correlation  between  the 
spFRET  and  the  open-closed  states  of  the  heterodimer 
gramicidin  strongly  suggest  the  existence  of  more  than  one 
conformational  state  at  each  open  and  closed  state  of  the 
gramicidin  ion  channels  in  the  lipid  bilayer.  These  hidden 
multiple  conformational  states  are  not  observable  from 
ensemble-averaged  measurements,  and  also  are  “silent”  and 
not  resolvable  in  a  conventional  patch  recording  analyse 


Homodimers  of  dye-labeled  gramicidin 
statistically  coincide  with  maximal  fluorescence 
self-quenching  at  an  open  state 

We  have  obtained  fluorescence  images  of  single  TMR- 
gramicidin  homodimers  at  bilayer  patches  simultaneously 
with  measurements  of  their  single-channel  electric  currents 
using  PCFM.  A  segment  of  a  single-channel  electric  current 
trajectory  (Fig.  5  a)  with  synchronous  fluorescence  images 
(Fig.  5  b)  illustrates  the  correlation  between  the  transitions  of 
the  channel  from  the  closed  to  the  open  state  and  fluo¬ 
rescence  hot  spots  (Fig.  5).  Single  hot  spots  were  analyzed 
for  width  (Fig.  6  a)  and  intensities  (Fig.  6,  b  and  c).  These 
mean  values  are  from  single  homodimers  of  die  dye-labeled 
gramicidin,  an  interpretation  based  on  the  coincidence  of  the 
recorded  single-channel  electric  currents  with  the  fluores¬ 
cence  images  and  on  the  various  control  experimental 
analyses  discussed  above. 

We  further  observed  that  the  width  of  the  hot  spot 
fluorescence  increased  with  fluorescence  intensity  (Fig.  6  a). 
Theoretically,  the  spread  functions  of  the  one-point-source 
and  two-point-source  images  are  Gaussian  (G(x,y))  =  A 
exP(  [  (x— Xq)2  +  (y-y0)2]/B }  and  a  convolution  of  two 
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FIGURE  5  PCFM  measurement  of  c 
single  TMR-grarmcidin  homodimer  :: 

ion  channels,  (a)  Upper  trace:  single- 
channel  electric  current  trajectory  re¬ 
corded  at  50-mV  polarization  potential 
from  a  single  TMR-gramicidin  channel 
at  a  patch-clamp  pipette  tip.  Lower 
trace:  time  flags  indicating  the  5  ms  of 
laser  excitation  and  CCD  camera  expo¬ 
sure  simultaneously  acquired  with  the 
single-channel  electric  recording.  “C” 
and  “0”  denote  the  electric  current  cor¬ 
responding  to  closed  and  open  states, 
respectively,  of  the  single  channel. 

(b)  Top  panel:  Consecutive  fluorescence  images  of  the  pipette  tip  taken  simultaneously  with  the  single-channel  electric  current  recording  shown  in  a. 

Bottom  panel:  plot  of  the  2-D  Gaussian  fitting  for  the  above  images.  Right  inset:  3-D  plot  of  the  first  image  (top)  and  3-D  plot  of  the  2-D  Gaussian  fitting 
(bottom)  with  subtracted  background  intensity  added  to  the  base  shown  in  gray. 
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Gaussian  functions  ( G  X  G),  respectively.  Assuming  that  the 
imaging  spatial  resolution  is  optical  diffraction-limited,  the 
FWHM  is  a  constant,  ~2£0'5,  for  G(x,  y)  at  different  signal 
intensity  and  the  FWHM  of  G  X  G  is  dependent  only  on  the 
two-point-source  separation.  The  error  bars  for  data  points  in 
Fig.  6  a  are  contributed  from  the  shot  noise,  finite  signal-to- 
noise  ratio,  digitized  photon  counting,  and  digitized  pixels. 
The  widths  of  the  hot  spots  were  essentially  diffraction- 
limited  but  wider  when  the  channels  were  fully  closed  (Figs. 
5  b  and  6  a).  This  observed  behavior  is  beyond  the  standard 
deviation  (Fig.  6  a).  These  observations  are  consistent  with 
the  expectation  that  the  channel  is  fully  open  when  the  dimer 
is  fully  associated,  leading  to  maximal  fluorescence  self¬ 
quenching  because  the  two  labeled  dye  molecules  have  the 
closest  distance  across  the  lipid  bilayer  at  the  patch  (Fig.  4  a). 
Rarely  (<5%  of  images)  were  two  separate  and  spatially 
resolved  imaging  spots  observed  when  a  single  channel  was 
closed,  i.e.,  a  dissociated  dimer.  In  most  cases,  the  two 
dissociated  monomers  do  not  separate  beyond  the  imaging 
spatial  resolution,  ~300  nm,  or  even  beyond  the  diffraction 
limit.  The  imaging  signal-to-noise  ratio  is  not  high  enough 
for  obtaining  a  clear  image  of  two  separate  spots,  as  we 
observe  that  the  widths  of  the  spots  are  larger  than  the 
diffraction  limit  at  higher  fluorescence  intensity  (Fig.  6  a). 

We  found  that  the  diffusional  motion  of  single  gramicidin 
channels  at  the  dPhPErDPhPC  bilayers  was  confined  with¬ 
in  the  imaging  spatial  resolution.  According  to  viscosity 
studies,  dPhPE  is  likely  to  be  in  a  gel-phase  at  room 
temperature  (Janko  and  Benz,  1977).  Based  on  our  control 
measurements  in  a  large,  supported  LB -membrane  bilayer, 
the  diffusion  coefficient  of  TMR-gramicidin  was  estimated 
to  be  <10“ 10  cm2/s  in  a  dPhPE  LB  layer,  yielding  ~0.2-jum 
average  displacement  in  1  s.  This  is  consistent  with  the  fact 
that  we  rarely  observed  diffusion  of  a  dissociated  dimer  into 
two  resolved  monomers.  Although  the  diffusion  coefficient 
of  the  TMR-gramicidin  can  be  as  high  as  10“8  cm2/s  in 
DPhPC  at  room  temperature  (Borisenko  et  al.,  2003),  the  4:1 
dPhPE:DPhPC  bilayers  may  have  small  domains  (<500  nm) 


of  the  DPhPC  formed  (Heller  et  al.,  1998),  confining  the 
diffusional  motions  of  the  TMR-gramicidin  molecules 
within  the  DPhPC.  Such  “ ‘corralled”  motion  may  be  ac¬ 
countable  for  the  rare  observation  of  an  optically  resolved 
separation  of  a  pair  of  TMR-gramicidin  monomers  dissoci¬ 
ated  from  a  dimer.  It  has  been  known  (Heller  et  al.,  1998) 
that  domains  of  DPhPC  can  form  in  DPhPE  bilayers  and  that 
gramicidin  channels  typically  prefer  to  stay  in  the  DPhPC 
regions. 

The  diffraction-limited  single  hot-spot  images  of  TMR- 
gramicidin  dimers  also  reflected  different  degrees  of 
fluorescence  self-quenching  (Zhuang  et  al.,  2000).  The 
self-quenching  critical  radius  of  the  rhodamine  6G  dyes  is 
~46  A  (Penzkofer  and  Lu,  1986).  The  fluorescence  self¬ 
quenching  between  the  two  dye  molecules  causes  an 
intensity  decrease  as  the  two  monomers  become  closer  and 
form  a  conductive  channel.  The  average  distance  between 
the  transition  dipoles  of  the  TMR  dyes  was  calculated  to  be 
~50  A,  and  the  transition  dipoles  were  nearly  aligned  in  the 
fully  dimerized  state  (Fig.  4  a). 

The  mechanism  of  fluorescence  self-quenching  is  still 
largely  unknown  and  system-dependent  (Zhuang  et  al., 
2000;  Penzkofer  and  Lu,  1986).  In  our  gramicidin  ion 
channel  system,  because  the  two  TMR  dye  molecules  are 
not  in  contact  to  form  a  dimer  or  excimer,  the  possible 
mechanism  of  the  observed  weak  self-quenching  effect  is 
likely  associated  with  reabsoiption,  and  secondary  fluores¬ 
cence  depends  on  the  fluctuating  orientation  of  the  two  TMR 
dye  molecules  across  the  lipid  bilayer. 

Fig.  6,  b  and  c>  shows  the  distributions  of  the  fluorescence 
imaging  intensity  in  correlation  with  the  open  and  closed 
states  of  the  gramicidin  channels,  measured  by  PCFM.  The 
intensity  of  the  hot  spots  was  statistically  higher  when  the 
channel  was  closed  (Fig.  6  b\  and  lower  when  the  channel 
was  open  (Fig.  6  c).  The  means  of  the  distribution  were 
found  to  be  67  ±  8  counts/ms  at  closed  states  and  52  ±  7 
counts/ms  at  open  states.  Furthermore,  the  difference  in 
fluorescence  intensities  of  the  single-channel  images  taken 
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FIGURE  6  (a)  Fluorescence  intensities  of  single  hot  spots  versus 
fluorescence  imaging  width  of  the  spots  as  derived  from  2-D  Gaussian  fit. 
{by  c)  Distributions  of  the  single-channel  fluorescence  intensity  for  closed 
and  open  states  of  single  channels  with  PDF  and  fitting. 


for  the  open  and  closed  states  was  also  found  to  be 
statistically  significant,  based  on  the  paired  r-test  (t  =  2.58, 
P  =  0.0297).  We  also  noted  that  the  broadness  of  the 
distribution  is  significant  and  beyond  the  measurement  noise 
and  that  both  quenched  and  unquenched  gramicidin  channel 
dimers  probably  existed  in  the  channel  closed  state.  Because 
the  efficiency  of  self-quenching  depends  on  the  distance 


and  the  relative  orientation  between  the  two  molecules,  the 
broader  distributions  of  the  fluorescence  intensity  at  the 
closed  (Fig.  6  b)  and  open  (Fig.  6  c )  states  also  suggests  that 
the  two  monomers  form  “closed”  or  “open”  dimers  of 
different  conformations  (Fig.  4).  This  result  is  consistent 
with  the  spFRET  imaging  result  (Fig.  3),  also  suggesting  that 
there  are  multiple  intermediate  conformation  states  corre¬ 
sponding  to  different  distances  between  the  homo-pair  of 
TMR  molecules  in  a  single  gramicidin  channel  based  on  the 
wide  distribution  of  the  fluorescence  self-quenching  effi¬ 
ciency  among  the  single  gramicidin  channels  in  the  lipid 
bilayer. 


Autocorrelation  analysis  of  single-channel 
electric  current  trajectories  characterizes  the 
dynamics  of  conformational  changes  of  the 
multiple  closed  and  open  states 

To  further  explore  the  possible  existence  of  the  multiple 
conformational  states  involved  in  gramicidin  ion  channel 
dynamics,  we  applied  autocorrelation  function  analysis 
(Zwanzig,  1990;  Neher  and  Stevens,  1977;  Labarca  et  al., 
1985;  Lu  et  al.,  1998)  of  the  single-channel  electric  current 
trajectories.  Detailed  single-channel  current  analyses  support 
our  attribution  that  multiple  conformational  states  are  in¬ 
volved  in  the  spatially  confined  gramicidin  dimer  channel 
dynamics.  Fig.  7  a  shows  a  typical  single  gramicidin  chan¬ 
nel  electric  current  trajectory.  The  intermediate  conductive 
states,  in  addition  to  the  fully  open  and  closed  states,  are 
evident  in  the  trajectory  and  are  further  revealed  in  an  elec¬ 
tric  current  amplitude  distribution  shown  in  Fig.  7  b.  There  is 
a  significant  population  beyond  the  measurement  noise 
and  within  the  time  resolution  in  the  distribution  between  the 
typical  open-closed  electric  current  amplitude  bimodal 
distribution  (Fig.  7  b).  The  autocorrelation  ftmction  deduced 
from  the  single-channel  electric  current  trajectory  (Fig.  7  a) 
shows  nonexponential  decay  (Fig.  7  c),  and  the  nonexpo¬ 
nential  autocorrelation  functions  were  commonly  observed 
for  both  dye-labeled  and  non-dye-labeled  gramicidin  chan¬ 
nels  (D.  J.  Panther,  G.  S.  Harms,  G.  Orr,  and  H.  P.  Lu, 
unpublished  results).  The  nonexponential  autocorrelation 
functions  suggest  that  the  transitions  between  multiple  open- 
closed  states  of  a  single  channel  either  occur  at  a  fluctuating 
rate  (Zwanzig,  1990;  Xie,  2002;  Yang  and  Xie,  2002;  Jung 
et  al.,  2002)  or  involve  a  complex  kinetic  mechanism 
associated  with  multiple  intermediate  states  (Schenter  et  al., 
1999;  Jung  et  al.,  2002;  Reilly  and  Skinner,  1993).  We  have 
analyzed  the  dynamics  using  an  autocorrelation  function  of 
sequential  open  dwell-time  (index  number,  m)  {ton(m)}, 
calculated  by  (Aron(0)A?on(m)),  where  A ton(m)  =  ton(m)  - 
(ton(ttO)  (Fig.  7  d).  The  single-exponential  decays  that  were 
fitted  to  the  autocorrelation  functions  of  ton(m)  (Fig.  7  d)  and 
toffcn)  (not  shown)  were  2.6(±1.0)  X  10"2  m~l  and  0.5  ± 
0.2  respectively.  The  open  and  closed  dwell-time 
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FIGURE  7  (a)  Single-channel  electric  current  trajectory,  { j*(t) } ,  recorded  at  100-mV  polarization  potential  from  a  single  gramicidin  channel  in  a  lipid  bilayer 
patch  (4: 1  dPhPE:DPhPC)  formed  at  a  patch-clamp  pipette  tip.  The  single-channel  electric  current  trajectory  indicates  that  there  are  only  single-channel  events 
(O’Connell  et  al.,  1990).  ( b )  Electric  current  amplitude  histogram  calculated  from  the  single-channel  trajectory  in  a.  Intermediate  electric  current  amplitude 
reflecting  the  inhomogeneous  conductivity  is  evident  although  the  channel  shows  essentially  two-state  kinetics.  “C”  and  “O”  denote  the  current 
corresponding  to  closed  and  open  states,  respectively,  (c)  Autocorrelation  function,  C(f)  =  (Az(O)Az(r)},  calculated  from  the  single-channel  trajectory  in  a . 
Biexponential  fitting  gives  kx  =  4.0  ±  1.0  s  ~\AX  =  0.4S,  and  k2  =  53  ±  5  s"\  A2  =  0.52;  using  Ax  exp  (-kit)  +  A2  exp(-*2r).  (d)  Autocorrelation  function  of 
“open”  dwell  times,  r{m)  =  (At(0)At(m)),  where  m  is  the  index  of  the  “open”  dwell-time  trajectory  for  the  single  TMR-gramicidin  ion  channel  recorded  in  a. 
The  single-exponential  decay  rate  constant  is  2.6  (±  1.0)  X  10-2  m”\  and  the  first  point  (not  shown)  of  the  autocorrelation  contains  the  uncorrelated  noise  and 
fast  fluctuation  beyond  the  instrument  resolution.  Error  bars  are  shown  for  the  first  ten  data  points.  The  calculated  autocorrelation  times  are  in  the  timescale  of 
seconds  and  subseconds  (c  and  d).  The  channel  open-closed  dwell  times  are  also  at  subsecond  timescale  (a).  Using  a  4-kHz  filter  instead  of  a  6-kHz  filter  does 
not  affect  the  results,  even  if  the  possible  aliasing  effect  exists  at  millisecond  timescale.  Inset:  The  distribution  of  “open”  dwell  time  deduced  from  the  same 
trajectory  in  a.  The  solid  line  shows  a  distribution  envelope  fitting  with  a  bi-Gaussian  function  with  a  relative  weight  of 0.808  and  0. 192  at  averaged  dwell  time 
of  1.13  X  10"1  and  8.67  X  10  3  s,  respectively.  Inhomogeneous  distribution  is  evident 


correlations  suggest  that  the  channel  open  and  closed  states 
are  dependent  within  a  short  period  of  “memory”  time  (Lu 
et  al.,  1998;  Zwanzig,  1990;  Yang  and  Cao,  2001;  Vlad  et  al., 
2002;  Jung  et  al.,  2002;  Xie,  2002;  Yang  and  Xie,  2002).  The 
memory  time  for  the  open  dwell-time  was  determined  to  be 
~8  s  (38  cycles),  but  only  ~0.4  s  (8  cycles)  for  the  closed 
dwell-time.  We  postulate  that  the  “memory”  effect  orig¬ 
inates  from  the  channel  open-closed  conformational  motions 
that  are  trapped  in  metastable  processes  under  different  rates 
of  toggling  among  different  conformational  states.  This 
dynamically  inhomogeneous  behavior,  also  defined  as 
dynamic  disorder  (Zwanzig,  1990),  has  been  observed  in 
other  ion  channel  systems  (Sakmann  and  Neher,  1995;  Neher 
and  Stevens,  1977).  Recently,  single-molecule  spectroscopy 
has  made  progress  in  understanding  this  dynamic  disorder 


behavior,  which  is  beyond  the  scope  of  the  conventional 
kinetics  (Lu  et  al.,  1998;  Zwanzig,  1990;  Yang  and  Cao, 
2001;  Vlad  et  al.,  2002;  Jung  et  al.,  2002;  Xie,  2002;  Yang 
and  Xie,  2002).  Typically,  the  mechanism  of  the  dynamic 
disorder  can  be  analyzed  by  the  “Agmon-Hopfield” 
diffusive  model  (Agmon  and  Hopfield,  1983;  Agmon, 
2000)  or  by  the  multiple-state  Markovian  models  (Schenter 
et  al.,  1999;  Sakmann  and  Neher,  1995;  Neher  and  Stevens, 
1977).  The  Fig.  7  d  inset  shows  the  distribution  of  the 
“open”  dwell-time.  At  least  two  subgroups  of  “open”  dwell 
times  contributed  to  this  broad  distribution.  We  postulate  that 
the  dynamic  disorder  of  this  gramicidin  system  in 
dPhPErDPhPC  bilayers  can  be  characterized  by  a  simple 
coupled  Markovian  kinetic  model,  for  example,  the  2  X  2 
model  (Schenter  et  al.,  1999),  which  describes  two  slowly 
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interconverting  open-closed  kinetic  processes  that  have 
different  kinetic  rates. 

The  open  and  closed  dwell  times  are  found  to  be 
dependent  over  periods  of  hundreds  of  milliseconds,  which 
suggests  a  trapped  or  not  fully  dissociated  subset  of  channel 
conformers.  We  note  that  the  gramicidin  channel  dynamics 
can  be  significantly  different  using  different  lipid  bilayers 
because  the  channel  open  and  closed  activities  are  associated 
with  the  gramicidin  diffusion  motion  in  a  particular  bilayer. 
Nevertheless,  the  gramicidin  ion  channel  dynamics  at  the 
dPhPE:DPhPC  bilayers  are  apparently  associated  with  more 
than  two  states:  1),  The  open  states  of  the  gramicidin  dimer 
involve  slow,  strongly  correlated  conformational  fluctua¬ 
tions  that  lead  to  dynamic  disorder  of  the  dissociation  rate 
(Lu  et  al.,  1998;  Zwanzig,  1990;  Yang  and  Cao,  2001;  Vlad 
et  al.,  2002).  2),  Gramicidin  dimerization  involves  multiple, 
weakly  correlated  or  noncorrelated  association  kinetics  at 
different  rates.  Although  the  association  kinetics  of  grami¬ 
cidin  have  been  studied  and  discussed  in  terms  of  fast 
fluctuations  (Sigworth  et  al.,  1987)  and  the  existence  of 
possible  “mini-channels”  (Busath  and  Szabo,  1981;  Sawyer 
et  al.,  1989),  our  observations  of  the  “memory”  effect  in  the 
gramicidin  channel  dynamics  in  dPhPErDPhPC  bilayers 
provide  a  new  insight  into  the  molecular  level.  However,  we 
are  not  sure,  at  this  stage,  if  the  complex  mechanism  we 
observe  is  generally  true  for  the  gramicidin  ion  channels  in 
other  lipid  bilayer  systems. 


A  multistate  model  for  gramicidin 
channel  dynamics 

Based  on  the  information  derived  from  the  PCFM  imaging, 
we  propose  a  working  model  of  multi-state  gramicidin 
channel  dynamics  at  dPhPE:DPhPC  bilayers  (Fig.  8).  The 
open  {left),  closed  {right),  and  intermediate  {middle)  states 
can  be  presumably  the  outcome  of  intra-  and  intermolecular 
conformational  changes,  geminate  recombination,  and  non¬ 
geminate  recombination.  Although  the  valine  (Townsley 
et  al.,  2001;  Sigworth  et  al.,  1987)  and  tryptophan  residues 
(Townsley  et  al.,  2001)  may  undergo  large  structural 
fluctuations  and  alter  channel  conductance  by  perturbing 
the  structure  of  the  dimeric  channel  (Koeppe  and  Andersen, 
1996;  Townsley  et  al.,  2001;  Koeppe  et  al.,  2000;  Salom 
et  al.,  1995;  Woolf  and  Roux,  1994),  we  postulate  that  the 
intermediate  conformers  result  predominately  from  the  six 
intermolecular  H-bond  fluctuations  that  disrupt  the  pore 
structure  of  the  gramicidin  dimer.  Conformational  changes  at 
the  N-terminal  linkage  of  the  dimer  may  significantly  dis¬ 
rupt  the  column  of  structured  water,  the  channel  cavity 
structure,  and  the  electrostatic  field,  thereby  perturbing 
channel  conductance  (Eisenberg,  1998).  Drastic  conductance 
modifications  were  reported  for  gramicidin  channels  with 
mutated  N-terminal  head-to-head  connections  (Sigworth 
et  al.,  1987;  Szabo  and  Uny,  1978),  indicating  that  the  six 
intermolecular  H-bonds  formed  at  the  contact  region  of  the 
dimer  stabilize  the  folly  conductive  gramicidin  channel 


FIGURE  8  Proposed  multistate  model  of  gramicidin  dynamics  indicating  a  third  “multistate”  of  nonconducting  or  partially  conducting  channels  with  larger 
separation  between  the  fluorophore  probes.  In  this  model,  the  third  “multistate”  represents  more  than  one  intermediate  conformation  state.  Alternative 
rotationally  distorted  conformations  are  possible  for  the  fully  conductive  state  of  the  gramicidin  channels  (Roux  and  Karplus,  1994). 
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(Szabo  and  Uny,  1978;  Koeppe  and  Andersen,  1996).  Partial 
disruption  of  these  H-bonds  results  in  incomplete  association 
or  dissociation  of  the  dimers  (Szabo  and  Uny,  1978).  Given 
that  each  H-bond  has  ~3-5  kcal/mol  dissociation  activation 
energy,  the  dissociation  of  these  bonds  plausibly  occurs  in 
the  millisecond  timescale  at  room  temperature.  These 
estimated  H-bond  dissociation  times  are  in  the  same 
timescale  as  the  gramicidin  dimer  changes  between  open 
and  closed  states.  The  conformational  changes  that  could 
originate  from  intramolecular  changes  or  lipid-channel 
interactions  are  expected  to  occur  at  a  faster  timescale  and, 
therefore,  should  not  result  in  the  observed  self-quenching 
and  FRET  effects  at  the  millisecond  timescale.  Therefore,  we 
hypothesize  that  the  conformational  changes  that  we 
measured  arise  from  the  N-terminal  H-bond  dynamic 
changes.  We  are  not  able  to  provide  direct  structural 
evidence  of  the  H-bond  fluctuation  using  our  current  PCFM. 
The  structural  and  dynamical  information  may  be  obtained 
by  integrating  a  Raman  spectroscopic  probe  with  the  PCFM, 
which  is  beyond  the  scope  of  this  article. 

Diffusion  of  individual  monomers  that  dissociate  from  the 
dimerized  state  (Tank  et  al.,  1982)  is  confined  within  the 
optical  diffraction-limited  spot  size,  based  on  our  experimen¬ 
tal  observation  using  wide-field  imaging.  This  observation 
has  been  further  confirmed  by  time-dependent  and  static 
anisotropy  analysis  of  gramicidin  C-TMR  in  the  LB  bilayer  at 
the  patch-clamp  tips.  The  detailed  results  and  discussion  will 
be  reported  elsewhere.  Briefly,  the  decay  time  of  the  ani¬ 
sotropy  is  ~5  ns,  and  the  static  anisotropy  is  0.25  ±  0.05 
(Harms  et  al.,  2003;  G.  S.  Harms,  G.  Orr,  and  H.  P.  Lu, 
unpublished  data).  Based  on  our  observation,  in  most  of  the 
events,  the  pairs  of  gramicidin  monomers  did  not  diffuse  away 
from  each  other  beyond  the  imaging  diffraction-limited  scale 
(^300  nm),  and  we  postulate  that  the  geminate  recombination 
may  be  associated  with  single-channel  conformational 
changes.  It  is  possible  that  the  diffusion  of  gramicidin 
monomers  was  confined  to  domains  (<300  nm)  likely 
associated  with  the  inhomogeneity  (Nielsen  and  Andersen, 
2000)  of  the  dPhPE:DPhPC  bilayers,  as  discussed  above. 

CONCLUSIONS 

We  have  developed  a  new  single-molecule  approach,  PCFM, 
combining  single-molecule  fluorescence  imaging  micros¬ 
copy  and  single-channel  patch  recording  to  simultaneously 
record  the  single-channel  trajectories  of  fluorescence  imag¬ 
ing  and  channel  electric  current  recording.  We  have 
demonstrated  the  application  of  PCFM  in  studying  the  ion 
channel  dynamics  of  single  gramicidin  dimers  in  a  lipid 
bilayer  at  the  patch-clamp  tip.  By  correlating  single-channel 
fluorescence  quenching  and  spFRET  imaging  with  single¬ 
channel  electric  current  recording  simultaneously  acquired, 
our  results  strongly  suggest  the  existence  of  multiple  con¬ 
formational  states  of  a  single  gramicidin  channel  associated 
with  the  channel  open-closed  activity  at  a  dPhPEiDPhPC 
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bilayer.  Although  the  inhomogeneous  dynamics  of  the 
conformational  change  were  probed  in  a  lipid  bilayer,  the 
exact  structures  of  the  conformational  states  require  further 
studies,  for  example,  by  correlated  Raman  spectroscopy.  The 
new  single-molecule  experimental  approach  for  studying  ion 
channel  conformational  dynamics  and  mechanisms  should 
be  applicable  for  direct  study  of  ion  channel  proteins  in 
living  cells  (Orr  et  al.,  2002). 
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